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Introduction 

1. The w-body problem is the name usually given to the 
problem of the motion of a system of many particles 
attracting each other according to Newton’s law of gravita¬ 
tion. This is the classical problem of mathematical nat¬ 
ural science, the significance of which goes far beyond the 
limits of its astronomical applications. 

The n-body problem has been the main topic of celestial 
mechanics from the time of its inception as a science. 
Now that many of the problems of celestial mechanics 
have become part of geophysics, the central position of 
the n-body problem has been further strengthened. 

The fundamental dynamical problem for a system of n 
gravitating bodies is the investigation and predetermina¬ 
tion of the changes in position and velocity that the 
particles undergo as the time varies. However, this 
is a very complex non-linear problem whose solution has 
not been possible under the present-day status of math¬ 
ematical analysis. 

It is natural that first to be developed were the prac¬ 
tical aspects of the n-body problem, arising in connection 
with the study of the motion of the planets, comets, and 
satellites. These problems were solved by the method 
of successive approximations. However, the mathematical 
difficulties entailed by these problems were so great, that 
the convergence of the approximations was not proved 
nor an estimate of the error obtained. Although the 
computational methods were not rigorously justified, their 
applicability, nevertheless, has been verified by practical 
usage which has shown that there exists an amazing agree¬ 
ment between theory and observation in most of the cases 
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considered. Despite their mathematical imperfection, the 
practical methods of celestial mechanics yield an accurate 
picture of the motion of the celestial bodies over a limited 
interval of time. 

2. Today, the hope for achieving some success in the 
study of the n-body problem is centered mainly on the 
use of qualitative methods of investigation. 

A qualitative investigation of the w-body problem 
usually consists in studying an incomplete set of parameters 
which, though not defining the state of the system uniquely, 
nevertheless, characterizes it in certain respects. Clearly, 
the problem of the predetermination of the values of an 
incomplete set of parameters cannot be well-posed, and 
we must confine ourselves to ascertaining the general laws 
by which these parameters change. This restricted for¬ 
mulation of the problem, characteristic of qualitative 
investigations, has the advantage that a knowledge of all 
the integrals of the motion is not needed for its solution, 
and significant results can be obtained by using only the 
known integrals, or by analyzing directly the equations 
of motion and their various analytical consequences. 

However, in choosing an incomplete set of parameters, 
it is important that we choose a representative set capable 
of reflecting the essential characteristics of the motion, 
and which leads to meaningful ideas about the dynamical 
system. To the set of such parameters belong, first of all, 
the mutual distances between the bodies of the system. 
It is therefore not accidental, that the search for the most 
general properties of the distances between the bodies as 
the time tends to infinity (positively or negatively) has 
occupied a prominent place in the qualitative investiga¬ 
tions of celestial mechanics. This problem is often called 
the “problem of final motions”. 
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The first general and very simple results on the final 
motion of it gravitating bodies were obtained by Jacobi 
in 1842. However, after Jacobi’s work, attempts were 
made to study only the special case of three bodies. A 
systematic investigation of the final motion of three 
gravitating bodies was begun about thirty years ago by 
J. Chazy. For this case, he attempted to classify the 
types of motions possible and to investigate their connection. 

The results of Chazy proved to be very interesting, 
although one of his most important ones, the impossibility 
of capture, is incorrect. After O.Yu. Schmidt showed, 
with the help of the techniques of numerical integration, 
that capture is possible, the need arose to carry on a more 
detailed study of the final motion in the three-body 
problem. It was necessary to obtain effective qualitative 
theorems that give sufficient conditions for the occurrence 
of specific types of motions in the form of restrictions on 
the initial data. A general approach to the solution of 
this problem was presented by the author in the years 
1948-1951, and he obtained the very first criteria for the 
occurrence of hyperbolic and hyperbolic-elliptic type mo¬ 
tions in the three-body problem. 

Research in this direction was continued by G.A. 
Merman, K.A. Sitnikov, V.F. Proskurin, N.G. Kochina, 
G. Ye. Khrapovitskaya, and O.A. Sizova. Thus, by 
today, a considerable amount of scientific work on this 
question has been accumulated (see the bibliography at 
the end of the book). 

We stated above that after Jacobi’s investigations, no 
attempts were made to obtain new results on the character 
of the final motion in the general case of an arbitrary 
number of bodies. The mathematical difficulties of this 
problem have apparently caused researchers to stay away 
from it. Nevertheless, we believe that progress in the 
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general case is possible, although it seems it may not be 
rapid. 

Below, we are going to consider some qualitative 
methods of analyzing the w-body problem, as well as some 
of the results obtained by means of these methods. Of 
course, the totality of these results does not, by far, give 
a complete qualitative picture of the motion of n grav¬ 
itating bodies, and its clarification will require still further 
work. 

3. Let us give a brief outline of the content of the book. 

In the first chapter, we give the equations and general 
integrals of the rc-body problem, and we study the simplest 
theorems on the final motion due to Jacobi. 

In the second chapter, we consider means of applying 
the method of dimensional analysis to the n-body problem. 
As far as we know, dimensional analysis has not been used 
before in the investigation of the final motion; though it 
does not yield definitive results in this area, it is very use¬ 
ful in carrying out a preliminary analysis of the problem. 

In the third chapter, we present our “method of con¬ 
tinuous induction/ 1 and we consider some applications of 
it in which the final motion in the w-body problem is 
analyzed. This method allows us to obtain effective 
qualitative results, namely, it allows us to formulate suffi¬ 
cient conditions for the occurrence of certain types of 
final motions in the form of conditions on the initial state 
of the dynamical system. 

The fourth chapter is devoted to the method of 
invariant measure. The application of this method to 
the many-body problem has necessitated working out a 
number of theorems on the measure theory of dynamical 
systems. These are presented at the beginning of the 
chapter. At the end of the chapter, we prove some very 
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general theorems on the motion of a system of gravitating 
bodies using the method of invariant measure. 

In the fifth chapter, an attempt is made to analyze 
some cases of the evolution of a system of n gravitating 
bodies on the basis of celestial mechanics. Here, we shall 
be concerned with the processes which are of cosmogonical 
interest, and which are accompanied by the conversion of 
mechanical energy into non-mechanical forms. 
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CHAPTER 1 


Equations and General Integrals 
of the n-Body Problem. 

Simplest Theorems Concerning the Final Motion. 

1. We shall consider a system of n particles Pi, P 2 , 

. . . , P w with masses m h m 2i . .. , m n attracting each other 
according to Newton's law of gravitation. We shall study 
the relative motion of these particles. At first, we consider 
the motion of Pi, P 2 , . . . , P* with respect to the center 
of mass of the system. 

Let us take a set of rectangular coordinates with origin 
fixed at the center of mass of the system and with coordinate 
axes having fixed directions. Let x i} y it Zi be the coordi¬ 
nates of the point P t *, r,- the distance of P,* from the origin, 
and tij the distance between Pi and Pj. We shall use a 
system of units for which the universal gravitational 
constant is 1. 

Let us set 


Jim ns 


( 1 . 1 ) 


where the summation is taken over all pairs of distinct 
bodies; then the equations of motion are given by the 
following 3 n differential equations of second order: 
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d*x t dU _ d*yt dU „ d*u dU 

m *SF=-Wi' m< rf<r == SST 

i 3 li tf • i •» M 


( 1 . 2 ) 


(In this chapter, we shall only consider the equations of 
motion and their integrals.) These equations may also be 
written as 6n equations of first order 


dxj 

fit 


— X 2 



dU 

dt 






dx i _ dU_ d v'i _ dU 

dt """ dxi * m i ffr • 








<= 1.2 .. 


(1.3) 


A first integral of these equations is the familiar law' 
of conservation of energy 


f '2 l m <(*? + y? + z‘ i 2 ) = U-t-H. (1.4) 

<=f 


where H is the constant of integration. 

Besides this, there are the six linear momentum in 
tegrals. In our choice of coordinates they are given by 



There are still the three area (moment of momentum) 
integrals which, in terms of the coordinates under con¬ 
sideration, are 
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s «i (v;—»!*;)= 

i=i 


S m «(yX— (L6) 

•=l 

2 m < (*X-*X)= c 3 . 

<=i 

Here Ci, £ 2 , and c* are constants of integration. 

These ten integrals exhaust all of the known essentially 
independent integrals of the given dynamical system. 

2. Another theoretically useful way of describing the 
relative motion of n bodies can be given by introducing 
Jacobi coordinates. We first change the notation of the 
preceding section into one that is better suited to the 
application of Jacobi coordinates. 

Let the particles in question now be denoted by P 0 , 
Pi,..., P»_i and their respective masses by Wo, mi, .. ., 
#V-1- 

Let £ 1 , rji , f i, be the coordinates of Pi with respect to 
axes with origin at the point P 0 , £ 2 , ^ 2 , £2 the coordinates 
of Pt with respect to axes with origin at the center of mass 
of the particles Po and Pi etc., and finally, let {*- 4 , 
^ 1 , 1 be the coordinates of the point P*— 1 with 

respect to axes with origin at the center of mass of the 
particles Po, Pi, ... f P rt -a- The coordinate axes in all 
of the systems are assumed to be parallel. Let p,- denote 
the distance of P< from the origin of the coordinate system 
in which the motion of this particle is being described. As 
before, we let r ii denote the distance between P,* and Pj, 
and we again use a system of units in which the universal 
gravitational constant is 1 . 
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The equations of motion in Jacobi coprdinates are then 
given by 

d*U dV d*r, d*U dU 

"rfTa —"517 ’ dt* — di )( ' ^ dt* ~ dCi * (1.7) 

* = 1, 2, .*i-l 

where /*i, n 2 , . . . , m *—1 are the reduced masses defined by 


”* 0 m T 

m o + m i * 


foo 4 w 1 ) ^2 
^0 4“ m i H" m 2 


u, _( m 0 + + . . . 4- ™*—o) 

" _1 m o + m \ -h • • • -+ m *_i 

Equations (1.7) may be expressed as 6(n —1) equations of 
the first order: 


d*j _ m _ • dCi _ 

dt dt ~ ^ rf/ — ^ 

*)<__ dU_ 

**' dt — dU ' (A< dt — </ij< ‘ dt~ (Kt * 

< = 1, 2, H—1 

In Jacobi coordinates, the energy integral becomes 

12 ‘‘.( s .' , +':’+C)-=^+«. ( 1 . 9 ) 


where H is the constant of integration. 
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The momentum integrals are satisfied identically in 
Jacobi coordinates, and the moment of momentum in¬ 
tegrals take on the form 


C* 7 !,)— c jt 

"|V,«-5,c;)= c ,. 


( 1 . 10 ) 


where Ci, c 2 , and c s are constants of integration. 


3, Finally, we give one more way of describing the 
relative motion of n gravitating bodies P 0 , Pi, ... , P„*-i 
using a coordinate system with origin located at one of 
the points, for instance, the point Po. 

The equations of motion in this case are given by 


• dt 2 


—m 0 



m J (E< — 6y) 




n —1 



dt 1 


~rn 0 


H—1 



N—1 



( 1 . 11 ) 


dKi 

dt* 






4+2 

r 0< It, 


i = i,: 


mj{ — Cy) 
r O 


n —1 


H-l 



where f < are the coordinates of P<. 

We shall also find it convenient to write these equations 
in vector form. Let i, j, k be unit vectors directed along 
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the £, i?, £ axes, respectively, and r 0 - the vector joining 
Pi to Pj . Then multiplying equations (1.11) respectively 
by 1, j, and k and adding, we obtain 


«-l n —1 



M, w—I 


We have just considered the motion of the system of 
particles relative to P 0 . However, we could consider the 
motion of all the particles relative to any particle P k . We 
may therefore write the following equations: 


**• 


* 1—1 




W—1 




m 4 


y=o 






1 , ...» *—1 

^“Op lp pi «p *•—1 


( 1 - 12 ) 


4. Let us now consider the moment of inertia of n 
gravitating particles: 





(1.13) 


Differentiating this twice with respect to the time t and 
taking into account that 




we obtain 
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2 + v<y'< + z / t '+*? 4- yf +O • 

(si 


By making use of (1.2) and (1.4), we can write this equa¬ 
tion as follows: 


II 



is *1 


d JL 

dxi 


+*%+*%)+*( 0 +B)- 


Now, the potential function U is a homogeneous function 
of the coordinates of degree —1; therefore 


2 f dV , dU , dU 
. (*<*+*■ «■+*■ 


iwm\ 


dii 


)=-« 


and hence, 


‘#=2(U + 2H). (1.14) 

This equation is called the Lagrange-Jacobi equation and 
plays an important role in celestial mechanics. It was 
obtained by Lagrange in the case of three bodies and was 
generalized to n bodies by Jacobi. 

5. We next derive an auxiliary equation which allows 
us to express the moment of inertia in terms of the dis¬ 
tances between the bodies. 

Consider the distances of the bodies Pi, Pa, . . • » P» 
from the center of mass of the system and the distances 
between the bodies. The square of the distance of P,- from 
the center of mass is 
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= (1.15) 


and the square of the distance between P { and P, is 

r 2 ij = (x i —x j ) 2 +(y i --y j ) 2 +(z i -~Zj) 2 . (1.16) 

We previously used to stand for the summation over 
all distinct pairs ij (i ^ j and ji not considered to be 
distinct from ij); we now let denote the summation 
over all possible pairs ij (any i is taken in combination 
with any j ). 

We then have the following obvious equations: 


V m,mj (x t — x ,)* — (*<—— 

*J *? 

= j 2 2 (x< _ — 


J=1 


n n » n n * 

- i2 2 miTnjX *+ 422 ~ 2 

<=i y=i <=u j =i <3 j=i 


=¥ 2 m<x *+t 2 m ' x * — 2 m<x< 2 

•=1 j—\ tel 


where 


M = 2 • 

tel 




However, since 
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and 

\ ** • 

Wj'Ej == 0 , 

by the linear momentum integrals, we find from the above 
equations that 

» 

M2 m i£ — (1.17a) 

isil -ij 

In an analogous way, we obtain the formulas 

« 

M 2 Ttity'i — 2 ((/< — yj)\ (1.17b) 

<=1 <> ' ' 

m 2 M j z? = 2m<'w > (3 < —Zy) 2 , (1.17c) 

<=i ij 

Adding (1.17a), (1.17b) and (1.17c) and taking into ac¬ 
count (1.15) and (1.16), we find 

.1/ 2 i m < r! = 2 m < m /*>- (1.18) 

<—i W 

Let us introduce the notation 



» 




Then by (1.13) and (1.18), we have 



Thus, the Lagrange-Jacobi equation is expressible in the 
form 
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i^ = 2(l/ + 2tf). (1.19) 


6. The energy integral and Lagrange-Jacobi equation 
permit us to establish certain properties of the distances 
r *j as time becomes negatively or positively infinite, 
without the use of additional analytical tools, for the 
general case of a system of n gravitating bodies. It is true 
that in this way we obtain only the simplest theorems on 
the final motion which already were given by Jacobi. 
However, the great importance of these theorems is that 
they reveal the proper way to formulate the problem of 
final motions, as well as indicate what general course a 
more detailed investigation should take. 

We first consider the properties of the distances r ij 
in a system of n gravitating bodies with a negative total 
energy. 

Theorem 1.1. If the constant H in the law of conservation 
of energy is negative , and if as t—y oo (t—y — oo J no collisions 
occur between the bodies , then for all t> 0 (/<0) the function 

r*(0-min{r^} 

ij 

possesses the following properties: 

(1) For all t>0{t<0) 




A/* 

2 | H\' 


(2) For any T>0 no matter how large and any 6 >0 no 
matter how small , there exists a time t’ with \t'\>T such 
that 
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<0< 


,V/2 


4 | U | — 2t » 


where H 

.U=2m ( . 


Let us prove the first statement. Since the kinetic 
energy of the system cannot be negative, it follows from 
the conservation law that 

U + H ^0 


or 



numj 

— 


>\B\- 


We can strengthen this inequality and write 

-i— |. (1.20) 

* W V 

However, it is not difficult to see that 
2 m { mj < — • 

<T 

Thus without weakening inequality (1.20), we may write 


and hence 


M* 

2r, (f) 


>1*1. 


Af2 

2|/f| * 


M0< 


( 1 . 21 ) 
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We now prove the second statement. We show that for 
any positive T> no matter how large, we can find a t r > T 
for which 

U+ 2Il> — €. 


Assume the contrary. Then for some value of T we have 
U+2II^ -e (1.22) 

for all values of t> T. 

Integrating the Lagrange-Jacobi equation (1.19) twice 
with respect to time from T up to t> T, we obtain 

t t 

n (/) = / 2 (T) 4- I v (T)(t-T) + 2M Jf (U 4- 2 H) dt dt . 

v r 

and hence by (1.22) 

/ 2 (/) < P (T) + P' (T) ( t-T ) — Mb (t - Tf. 

We see from this that —a> as oo. Therefore 

there must exist a time r>0 such that l 2 (r) = 0, i.e., a 
collision between the bodies occurs. This contradicts the 
condition of the theorem. 

Thus, there exists a t f > T for which 

U(t')+2H >—€, 

U(t , )>2\H\-—e. 

Strengthening this inequality, we have 
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or 


%(/)< 


M 2 

41 jFZ | — 2 € ' 


and so for 0 the theorem is proved. For t^ 0, the 
theorem is proved in a similar way. 

The fact that when //CO, some of the distances can 
increase indefinitely with the time, does not contradict the 
theorem just proved. However, it follows from the theorem 
that this cannot happen to all the distances. Thus when¬ 
ever H< 0, it is impossible for a system of gravitating 
bodies to be completely scattered either as tr+ oo or as 
t — y — oo. 

We now consider a system of gravitating bodies for 
which the total energy is positive, i.e., //>0. 

Theorem 1.2. (Jacobi’s Theorem). If the constant H in 
the energy integral is positive , and if as t-+- oo (/->— no 
collisions occur between the bodies , then at least one of the 
distances Yij increases indefinitely . 

Proof: We start with the Lagrange-Jacobi equation 
in the form (1.19) and integrate it twice with respect to 
time from / = 0 to £>0. We thus obtain 

t t 

l*(t) — /» 4- Iff t- 2M Jf (U -A- 2H) dt dt, 

n o 


where II and Jo' are the respective values of P and I v 
at /=0. Now, since U>0, it follows that 

t t 

l'(t)>P, + P’t + AM\\Hdtdt, 

• 0 


or 


I*(t)>Il + P 0 t + 2MHt'. 
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When H> 0, we find from this inequality that 
P (/)-*» as £-»-+ 03 (<->-— 00 ), 

This is only possible if at least one of the distances r ti 
increases indefinitely as <->-4- 00 (<->-— °°). The theorem 
is proved. 

This theorem is not incompatible with the fact that, 
when H> 0, instances are possible in which all of the dis¬ 
tances between the bodies increase indefinitely as <->- 00 , 
or as <->-— ®, or else in both cases. However, the 
theorem is also not incompatible with the fact that, when 
H> 0, there may be subsystems in which the distances 
between the bodies remain bounded at <->-», or as <->•— 00 , 
or as t approaches both db 00 . Finally, the fact that there 
may be complete scattering of the system as the time 
increases in one direction coupled with incomplete scatter¬ 
ing as the time increases in the opposite direction is not 
incompatible with the theorem. 

In conclusion, by comparing Theorems 1.1 and 1.2, 
we can state that, when H< 0, the energy in the system is 
insufficient for it to be completely scattered, and when 
ff>0, the energy is so high, that at least a partial scatter¬ 
ing of the system is inevitable. 



CHAPTER 2 

Method of Dimensional Analysis 

1. The method of dimensional analysis, which has 
lately been applied rather widely in diversified problems 
of mechanics and physics, can also be used in the qualitative 
analysis of the motion of systems of gravitating bodies. 
However, it must be borne in mind that dimensional 
analysis is incapable of yielding complete enough results 
for the n-body problem, though in certain instances it 
enables important conclusions to be obtained. 

Our purpose is to consider means of applying the meth¬ 
ods of dimensional analysis to find criteria for the occur¬ 
rence of various types of motions in a system of n grav¬ 
itating bodies. In those instances where dimensional 
analysis enables us to obtain such a criterion, we find that 
it is expressible exactly to within an unknown dimension¬ 
less factor. Thus, the methods of dimensional analysis 
do not lead to a complete criterion, but only give an 
indication of its analytical form. Nevertheless, such re¬ 
stricted results are of importance since they yield pre¬ 
liminary knowledge and give us a correct orientation at 
the initial stages of our investigation. 

2. That dimensional analysis is capable of indicating 
the analytical form of the dependence of two dimensional 
quantities springs from the fact that objective physical 
relationships cannot depend on the choice of the system 
of units. However, such relationships between quantities 
are not arbitrary and are determined by the dimensions 
of these quantities. 

Dimensional analysis is usually applied to some rela¬ 
tion expressible in the form of an equation or identity. 

15 
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What is somewhat new in the application of the method 
of dimensional analysis to the qualitative study of the 
w-body problem is the necessity for considering inequalities. 
For our further purposes, the following reasoning will 
suffice. 

(1) Let SfK be a set of systems of basic independent units 


Ql) $ 2 > •••»<?* 

such that a transformation to any other system of units 
(7t> • • • » Qm 

also belonging to SfR, is effected by means of the formulas 

= ( 2 . 1 ) 
* - 1 , 2 . n 

where a,* are arbitrary positive numbers. 

(2) There exist arbitrary positive quantities 

Qoi Qh (?2, ..., (?«, 

which are homogeneous in the sense that the values of these 
quantities : 

V* Q'b Q'h . • • » Qn 

after transforming the basic units according to (2.1), are 
related to their previous values by the equations 

Q'j = afii a?i2... a&nQj. (2.2) 

The exponents a& are the powers of the dimensions; if 
they are all zero, then the Qj are dimensionless. It is clear 
that the product of homogeneous quantities is again 
homogeneous. 
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(3) The quantity Q 0 is related to Q h Q 2 , . . . , Q n by 
the inequality 

Qo^fiQu <? 2 . Qn\ (2.3) 

(4) The relation , Q n ) does not depend 

on the units of measurement in the sense it is satisfied after 
every transformation of the basic units by means of (2.1). 

(5) All the basic units enter into the expressions for 
the quantities Qo , (?i, . . . , Q„. 

Theorem 2.1. If conditions (1) through (5) are satisfied, 
then we have 


Q»^CQXiQ* 2 * (2.4) 


where the dimensionless constant C=/(l,l> . . . , 1), awd the 
exponents Xj are uniquely determined by condition (4). 

Proof: We take the quantities Q lt and 

transform to new units; we choose the scales ai,a,2, . • a* 
so that Qi, Qi, - - ■ , Qn each has the value 1 in the new 
system of units. Imposing this condition on (2.2), we 
obtain 


afji a\n • • • = 


1 

Q/ 


(2.5) 


Solving these equations, we find that the a j are given by 


a-’ 


1 

. . . Q*in ’ 


( 2 . 6 ) 


where bjk are constants whose values we shall not deter¬ 
mine. (We might suggest that the reader solve equations 
(2.5) for the simplest case of n=2 if he is interested in the 
details of the derivation of (2.6)). 

Let us now write (2.2) for the particular case where 

J= °' 
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Substituting in this the value of from (2.6), we obtain 

< 2 - 7 > 

whw-e for the moment we do not determine the exponents 

■*li *4» . . • , X n . 

Inequality (2.3) does not depend on the units nf 


Setting 


i.i). 

C-fi 1 . 1 . 1 ) 


"w inequaKy ,he «*—-'» 


and hence 


- -^° < r. 


toXdiXn ‘f4Wh thUS ^ ° btained ‘ Now -^cording 
■HmmienUy “ e<Iullit >' is 

obtain 

The theorem is proved. 

denote the distance from ?, ,, , nd der '^ 
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of this quantity with respect to time. We shall call the 
latter the radial velocity. 

Let us introduce the notation 

p(/) = min [r {j ( t)}, 
ij 

<r(t)= min {r' fi (t)}. 

*j 

We shall assume that at the initial time J = 0 the condition 

<r(0)>0. (2.8) 

is satisfied. 

Under this assumption, there are two possibilities: 

(a) it is possible to give a positive number r such that, 
for all i in the half-open interval [0,r), the distances 
between the bodies increase and from the moment of time 
r at least one of the distances ceases to increase; 

(b) or if cr(0) is so large, and therefore also the initial 

values rij CO) of all the radial velocities, that the mutual 
attraction of the bodies of the system is not great enough 
to contain the system, then all the r^(/) will approach 
infinity as /--> . 

However, for exactly what value of <r (0) will this 
phenomenon occur? This, of course, will depend on the 
initial state of the system, on the parameters of the system 
and on the parameters which characterize the phenomenon 
of gravitation in general. 

Let us try to clarify exactly which parameters are of 
great importance in the problem under consideration. 

From general physical notions about the force of 
gravity, it follows that the smaller the value of p(0), i.e., 
the closer the bodies are at the initial time, the greater 
must be the value of cr(0); therefore, there is no doubt 
that one of the parameters is p(0). 
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The masses mo,mi, ... , m w __i are also of considerable 
importance. For, the larger the masses of Pq, Pi, • • •»P«— 1 » 
the stronger will be their mutual attraction and, with all 
other conditions being equal, the greater must be the value 
of a ( 0 ). Thus, taking into account the masses will yield n 
more parameters. However, we can avoid the large number 
of parameters entailed by the masses if we remember 
that we are going to obtain sufficient conditions in the 
form of certain inequalities. Therefore for the purposes 
of simplifying the problem, we can use an overestimate of 
the role of the masses. As a parameter characterizing 
mass effects, we take 

M = mi, 

which is some function of the masses having the dimension 
of mass. 

Finally, gravity is of essential importance, and there¬ 
fore the gravitational constant 7 must also be one of the 
parameters of our problem. 

The parameters 7 , ju, p( 0 ) and <r( 0 ) constitute a 
complete set of parameters for the phenomenon under con¬ 
sideration. In fact, for given values of p( 0 ) and p, if 
tf( 0 ) is sufficiently great, then the bodies of the system 
will all be scattered, regardless of what values the other 
unaccounted for parameters may have. Thus we con¬ 
clude that the required criterion must have the following 
general form: 

<r( 0 )>/( 7 , M , p(0)), (2.9) 

where the function / is momentarily unknown. We must 
now determine the functional form of /( 7 ,p,p( 0 )). 

Consider arbitrary but fixed values of the parameters 
7 , p, and p( 0 ); let<r* ( 0 ) denote the greatest lower bound 
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of the set of values of a (0) for which the system of bodies 
is completely scattered. This lower bound obviously exists 
because of inequality (2.8). It is further obvious that 
the inequality 

<r*(0)^/( 7 , M ,p (0)) (2.10) 

is satisfied for all values of y, p, and p(0). 

Let us take the system of units ordinarily used in 
problems in mechanics in which the basic units are time T\ 
length L, and mass M. The parameters in which we are 
interested have the following dimensions: 


[P(0)]=L, 

[cr*(0)]-Z,r-L 

[y]= L 3 T~~ 2 M~~ l . 


( 2 . 11 ) 


We thus have four parameters, three basic units, and 
the inequality (2.10), i.e., our conditions correspond to 
those of Theorem 2.1, and to apply it, we must set n —3. 
By this theorem, we then have 

<r„c(0)^C 7 W(0)*, (2.12) 

where the exponents y, and z are for the moment un¬ 
known. To determine them, we equate the dimensions 
of the left-hand and right-hand sides of this inequality. 
This yields 


LT -1 = 


or 


LT~~ X =L 3x+t T^ 2x M y ^ x . 
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Equating exponents of like dimensions, we obtain the 
equations 

2x=l, 

y—x=0, 

3x-\-z= 1, 

whose solution is 

• E =y> y = - j» l — —y. 

Substituting these values of x, y, and z in inequality 
(2.12), we obtain 

(2-13) 


However, if our system of units is the one usually used in 
theoretical investigations and in which y=l, then the 
criterion becomes 




Our above result implies that 

/(T, |*, p(0)) = c )/3^, 

and thus (2.9) becomes 

’<°» c l / rfe (7=D- 


(2.13a) 


(2.14) 
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The value of C and the form of the function expressing 
fx in terms of mo,Wi, .. . , and w n _ t cannot be determined 
by means of dimensional analysis. Thus dimensional an¬ 
alysis does not yield definitive sufficient conditons for the 
occurrence of the type of motions of n gravitating bodies 
in which we are interested; it merely determines the 
analytical form of the criterion. Later on, we shall see 
that the results obtained are confirmed by the direct 
analysis of the equations of motion. 

In conclusion, let us examine the methodological as¬ 
pects of our above investigations. 

We obtained inequality (2.9) on the basis of general 
physical considerations. This part of the investigation 
does not concern dimensional analysis. Figuratively 
speaking, inequality (2.9) is not within the 4 'competence' r 
of dimensional analysis, and it cannot be “responsible" 
for this inequality. What merely follows from dimen¬ 
sional analysis is that if (2.9) is in principle correct, then 
it has the form (2.14). The question as to whether or 
not inequality (2.9) is the correct form of a sufficient 
condition for the occurrence of different types of motion 
cannot in principle be answered by dimensional analysis. 

In this connection, it is interesting to note the follow¬ 
ing. Our reasoning cannot be applied to the case where 
<r(0)<0. In fact, if we wrote the inequality (2.9) in 
the form 

I(0)|>/(7* M* P(0)), (2.15) 

we could, at first glance, repeat our argument and obtain 
a similar but no more general result. However, in reality, 
this is not so. For, in assuming that <r (0)<0, we would 
not have the right to suppose that the problem is solved 
by inequalities (2.9) and (2.14). When <r(0)<0, it is 
possible for the gravitating bodies to approach one another 
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very closely, and the resulting gravitational interactions 
could then lead to such an exchange of energy between the 
bodies that a portion of them could go off to infinity with 
increased velocity, whereas the energy of the remaining: 
bodies would be insufficient for them to be scattered. In 
this case, it is not sufficient to consider merely the 
parameters <r( 0 ), p( 0 ), ju and y in formulating conditions 
for the scattering of gravitating bodies, and we require 
a more detailed set of parameters that characterizes more 
fully the initial state of the system of gravitating bodies. 
Therefore, when <r( 0 )< 0 , the relation (2.14) must either 
be replaced by one containing a fuller set of parameters, 
or it must be supplemented by relations containing other 
necessary parameters, or else our considerations must he 
limited to such a set of initial states for which the un¬ 
desirable close approaches cannot occur. 

Let Vij be the magnitude of the velocity in the relative 
motion of P, and P jy and let 

z</)=min{y, v (0}. 

It is clear that if motions of gravitating bodies are to 
occur for which p (/)-> «> as /->■ oo , it is not at all necessary 
for the radial velocities to be sufficiently great. Such 
motions can occur if the absolute velocities in the relative 
motion of the bodies are sufficiently great. Nevertheless* 
we cannot substitute v(0) for a( 0 ) in inequality ( 2 . 9 )* 
even though the dimensions of these quantities are identi¬ 
cal. This is due to the fact i>( 0 ) is always positive, whereas 
it is possible for z>( 0 ) = | a (0)| when a (0) < 0. But we saw 
above that a criterion of the form ( 2 . 9 ) does not apply to 
this particular case. Therefore, any criterion for the 
occurrence of motions in which (/)-*■ oo as /-> co cannot 
contain merely p( 0 ), p( 0 ), p, and 7 , but must contain 
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still other parameters that, to one extent or another, 
explicitly or implicitly account for the direction of the 
velocity. As we shall see later on, this may be achieved 
by using fl(0) and <r (0) together since their difference, to 
a certain extent, characterizes the direction of the velocity. 
In this way, one may obtain a criterion which is applicable 
for a sufficiently wide range of values of cr(0)<0. 

We have limited ourselves to a comparatively simple 
application of the method of dimensional analysis. How¬ 
ever, this application makes rather well-apparent its 
characteristic features and the possibilities of its use. 

In summing up, we must admit that the method of 
dimensional analysis requires great care and thought, as 
well as a deep and penetrating intuition into the essence 
of a problem. Nevertheless, dimensional analysis can 
be earnestly recommended if it is simply regarded as an 
exploratory or research means in theoretical considerations. 




CHAPTER. 3 

Method of Continuous Induction 

1* We shall consider a set 3TC whose elements are time- 
dependent processes or mechanical or physical systems S 
whose states vary with the time t. 

Let r be some property of 5. In what follows, we 
shall make use of the following definitions. 

(1) The property T is said to be inductive on the set 911 
if for every S69H we have the following: from the assump¬ 
tion that the property V is satisfied for all values of t 
such that 

0 <Kr , 

it follows that it is satisfied for t=r no matter what r is. 

(2) The property T is said to be continuous if from the 
condition that it is satisfied at a certain moment of time 
Ji, it follows that it is satisfied in a certain neighborhood 
of h f i.e., it is possible to find a number 5 >0 such that Y 
is satisfied for all t such that 

I t—h\<8. 

Theorem 3.1. (On Continuous Induction). Let the 
property Y be continuous and inductive on the set 9TI. If Y 
is satisfied at the initial time /= 0, then it is satisfied for 
all values of />0. 

Proof: We call the time t regular if for all /' such that 

0 ^'^, 

the system S possesses the property T. Let Q denote 
the set of all regular times. From the continuity of V 
and the assumption that V is satisfied at /=0 T it follows 
that the set Q is not empty, 
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We prove the theorem by contradiction, and we assume 
that T is not satisfied for all values of Therefore, 

Q is bounded from above; let t = sup Q. Since t is the 
least upper bound, it follows that t is satisfied for all 
values of t in the interval 0 ^t<r. But then since T 
is inductive, it is satisfied at time r. Furthermore, the 
property T is continuous. Therefore, we can find a 
5 >0 such that T is satisfied for all / for which 

\t—r\< 5, 

But then V must be satisfied for values of t for which 



Hence, it follows that the time 

7" i = T —j-> T 

2 

is regular. However, this contradicts the fact that r is 
the least upper bound of the set Q . The theorem is proved. 

2. To apply the Theorem on Continuous Induction, 
we must have criteria for when some particular property 
is actually inductive. We shall now show that such 
criteria can be given in certain sufficiently general cases. 
Let 

COi, 0>2, • • . , W, 

be some set of parameters of the system 5 depending on 
the time which, while not completely determining the 
state of the system, characterizes it in certain respects. 
We denote the set of functions 


G>1, <02, • • W, 

by {co}. Moreover,let the system depend on ju parameters 
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which are independent of the time, and which we denote 
collectively by {a}. 

Let u4({a},{w}) denote an operator which, according 
to certain rules, sets in correspondence to each set of 
values of the parameters Oi, 02 ,. . and to each set of 
v functions coi, w 2 , * a continuous function ^(t) 
defined for all t ^ 0. 

We first consider the case where cui, w 2 , are 

unknown functions, and our knowledge of the system S 
is restricted to the information that these functions and 
the parameters ai, as, . . a M satisfy relations of the form 

«,(*)£ii,({«} f M). (3-U 

, 2 .*- 

Here, are operators which are defined by means of 
their analytical expressions. 

Let 

/i(0.A(*)t .--./.(O 

be suitably-chosen continuous functions of t defined for 
all O0 which we shall call comparison functions. The 
property T will be that the following inequalities are 
satisfied: 

» - 1 , 2 , ..«,<* 

Theorem 3.2 (First Theorem on the Inductiveness of 
Inequalities). // a system S satisfies the operator inequalities 
(3.1), and if for any 0, from the assumption that the 
inequalities 

»<«>/<( 0 
i - 1 , 2, .... »» 


(3.2) 
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We prove the theorem by contradiction, and we assume 
that T is not satisfied for all values of 0. Therefore f 
Q is bounded from above; let r = sup Q. Since r is the 
least upper bound, it follows that r is satisfied for all 
values of t in the interval 0^2<r. But then since T 
is inductive, it is satisfied at time t. Furthermore, the 
property T is continuous. Therefore, we can find a 
5>0 such that T is satisfied for all t for which 

|f—r|< 8 , 

But then T must be satisfied for values of t for which 

0^ T + -. 

2 

Hence, it follows that the time 

Tl=T-\ ->T 

2 

is regular. However, this contradicts the fact that r is 
the least upper bound of the set Q. The theorem is proved. 

2. To apply the Theorem on Continuous Induction, 
we must have criteria for when some particular property 
is actually inductive. We shall now show that such 
criteria can be given in certain sufficiently general cases. 
Let 

to%, . . ., oi v 

be some set of parameters of the system »S depending on 
the time which, while not completely determining the 
state of the system, characterizes it in certain respects. 
We denote the set of functions 

W 2 , . . ., CO, 

by {o>}. Moreover, let the system depend on parameters 
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&l> Os, • • ■ » Oji, 

which are independent of the time, and which we denote 
collectively by {a}. 

Let -4({a},{<*>}) denote an operator which, according 
to certain rules, sets in correspondence to each set of 
values of the parameters a%, <k>- . a* and to each set of 
v functions wi, c*> 2 , a continuous function ^(/) 

defined for all 0. 

We first consider the case where wi, o> 2 , ..•»«* are 
unknown functions, and our knowledge, of the system S 
is restricted to the information that these functions and 
the parameters a u 02 , . . satisfy relations of the form 

«<(0 ZA t <{«},{«}). (3-D 

»-l, 2-, * 

Here, A , are operators which are defined by means of 
their analytical expressions. 

Let 

. /,(/) 

be suitably-chosen continuous functions of t defined for 
all t^O which we shall call comparison functions. The 
property T will be that the following inequalities are 
satisfied: 

«, (0 >/, (<)• 

. ..* 

Theorem 3.2 (First Theorem on the Inductiveness of 
Inequalities). If a system S satisfies the operator inequalities 
(3.1), and if for any r^O, from the assumption that the 
inequalities 

«,(/)>/.(/) 
i -1 . 2. » 


(3.2) 
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are satisfied in the half-closed interval 0 ^ /< r, it follows 
that the inequalities 

A { ({a}), {»})}A { Ua}, {/}), (3.3) 

are satisfied , then (3.2) is inductive on the set 9TI of systems 
S for which the inequalities 

A-(W, {/})-/,(*) >0. (3.4) 

i-1,2. v 

are satisfied for all t^O. 

Proof: Let us write inequality (3.1) in the form 


(3.5) 

• “ l , 2. v 

Choosing an arbitrary r >0, we assume that 

«.(<)>/,(0 for <K/<t. (3.6) 

* -1, 2.v 

Then considering (3.5) for values of t in the interval 
0^<r and using (3.3), which holds because of (3.6), we 
can write 


({a}, {/}.)-/,•(<)• 

>-l, 2.. 

We now let t -> r taking into account (3.4) and the fact 
that /,•(<), and A ,• {{a ), {«}) are continuous func¬ 
tions of /, and we obtain in the limit 

«<(t)-/,-( r)^,.({a}, {/},)—/ < ( T )>0, 

*-l, 2, , v 


or 
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«<(T)>/,(r). (3.7) 

* -1 , 2 .* 

We have thus shown that the inequality (3.7) follows 
from (3.6), and this proves the theorem. 

3* We now consider a more complicated case of the 
Theorem on the Inductiveness of Inequalities which we 
shall prove in two parallel versions. 

Suppose we have two sets of parameters of the system 
S which depend on time: 


COl, CU2, ...» c Oy 

Si, 8 2» • • ., 8 V 

which, though not defining the state of the system com¬ 
pletely, characterize it in certain respects. As before, 
suppose the state of the system also depends on /i para¬ 
meters 4i, 02 , . . o P which are independent of time. We 
consider the case in which o>i, w 2 , ..c*>„ and 8 1 , $ 2 ,. . 8 y 
are unknown functions, and our knowledge of the system 
5 is restricted to the fact that these functions and the 
parameters Oi, a*, . . ., a M satisfy relations of the form 

{«}) 

M)). (3,8) 

i-l. 2, ...» * 

Here, A ,• are operators defined by their analytic expressions. 
Let 


be suitably chosen continuous functions of the time t 
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which are defined on the range 0, and which we call 
comparison functions. We shall assume that they pos¬ 
sess the following property: From the assumption that 
the inequalities 

5 « 00 < fi (0 

(3.9) 

<= 1 , 1 , ...» * 


are satisfied in the interval 0^/<r, it follows that the 
inequalities 


(0 >/<(*)• 

<= 1 , 2 .* 


(3.10) 


are satisfied for all values of t. 

Although we do not know the functions c<>,•(/) and 
*<W» we can, however, often establish such relations 
from the physical or mechanical meaning of these func¬ 
tions or from other considerations not requiring the know¬ 
ledge of the analytical form of these functions. 

We shall take Y to be the property that the inequali¬ 
ties (3.9) are satisfied. 

Theorem 3.3 (Second Theorem on the Inductiveness 
of Inequalities). If the system S satisfies the operator 
inequalities (3.8), and if for any r >0, from the assumption 
that the inequalities (3.10) are satisfied in the interval 
0 ^t<r it follows that the inequalities 


MW, WK^Ka), {/}) 

{*))> M{a), {/})), 
<=!.*.. 


(3.11) 
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are satisfied , then property (3.9) (and (3.10) which follows 
from it) is inductive on the set OH of systems S for which the 
inequalities 


A (W. {/}) — ?< (0<0 

U«a>,{/»-<!P,«)>0). 

<= 1 , 2 , 


are satisfied for all t ^ 0. 

Proof: Let us write inequality (3.8) in the form 

5 <(0— ?<(*)< M) — ?/(*) 

(s <(t)-?<(0>A i ((a}, {<*})-?,(*)) ( 3 - 13 ) 

*=1,2. 


Choosing an arbitrary r >0, we assume that 


2,(0 <?,( 0 . 

0</<t 

(M0>?< (0J- 


(3.14) 


Then for these same values of /, the inequalities 


«*('»/«(*)■ 
* = 1,2 . 


(3.15) 


will be satisfied. Considering (3.13) for values of t in the 
interval 0^t<r and making use of (3.11), which hold 
because of (3.15), we can write 
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Wenowlet t->T taking into account (3.12) and the fact 
that «,(<), 8 { (t), <p { (t), and -4<({a},{»}) are con¬ 

tinuous functions of t, and we obtain in the limit 


or 


5 <C0— ?<(*)<'UM. (/)) — 'fcC'XO 

(*< (*) - 9, (t) > A t «d), {/}) - * (t) > 0). 

•*» 1,2 . 


W < ?i (x) 

(s 4 (t)><p 4 (x)), 
#= 1,2 .. 


(3.16) 


From this, in turn, it follows that 

“><(*) >/<(*)• 

# = 1, 2, • ■ •! * 

We have thus shown that inequality (3.16) follows from 
(3.14) with which the theorem is proved. 

4* We now consider some applications of the Theorems 
on Continuous Induction and the Inductiveness of In¬ 
equalities to the qualitative analysis of gravitating bodies. 

We begin with the first Theorem on the Inductiveness 
of Inequalities. 

Let us take a system of n gravitating bodies 
Po, Pit . . P with masses m 0t m h . . i, and let 




rrumj 
m < + m / 


We shall make use of Jacobi coordinates. To begin with, 
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we consider the motion of Pi relative to Po. Differen¬ 
tiating the identity 


-2 — S2 ! r 2 , *2 

r Ql - "l I T| 1 + * 1 * 


twice with respect to t , we obtain 


r oi r oi + r oi “ + Wi + Vi + + \ + C’ 


*>2 


or 


t" i ’ll i r H l 

oi + C i + 


+ 


m ,2 t2 „»* 

+ T *1 +*•! “~ r 0l 


r 01 


Taking into account that 


.»2 , '2 . „4 io 

*1 +^l + t\ — r 01 
r 01 


>0, 


we find that 


r" "> it E* 1 31-11*4 ^ - T 


r oi 


r 01 


Using the first of equations (1.7), i.e., for 1, we obtain 


'oi 



r oi 


*U 


'‘oi 


Sl . 

roi 
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Now, let l denote the direction of the line joining P 0 to 
Pi. Then the last inequality may be written in the 
following form 

„ ^ 1 AU 

r oi^u 1ll oi * (3.17) 


Differentiating the potential function U in the direction 
l using (1.1), we have 


oU 

or 


—y 


v 


tmm j 


dry 

~dT 


From simple geometrical considerations, it follows that 



and, therefore, 



Inequalities (3.17) and (3.18) imply that 


(3.18) 


r" ^ 1 V m,m * 

ij o 


We have obtained this inequality by considering how 

•a m °. VeS relative to P °- However, we could have con¬ 
sidered the motion of any point P s relative to any other 
point P*. Thus, we have 


%*,>- — v 

* “M — 


rrtimj 

~ 


Finally, integrating this inequality twice with respect 
to time from 1-0 to *>0, we obtain the following system 

of integral inequalities: 
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r* (0 > ^ (0) + 4 (0)1t ■- ~ v j f 22L,dt dt. (3.19) 

-V n - ii 


ij 0 n 


We now use these inequalities to find sufficient con¬ 
ditions which ensure that the distance of each body from 
the remaining ones will become infinite as /->«, in other 
words t that the system is completely scattered. Such 
will be the case, at any rate, if the inequalities 


ruV)> 


r».( 0) 

~"2 



t 


are satisfied for all t ^ 0. 

To find the conditions, it is most convenient to use 
the first Theorem on the Inductiveness of Inequalities 
(Theorem 3.2). 

For the functions «,•(<), we take the distances r*,(0 
and for the parameters a it the initial values of these dis¬ 
tances and their derivatives, r*,(0) and r' fc ,(0), and the 
masses m 0 , mu .... wi»—i- 

In accordance with inequality (*), it is natural to 
choose the comparison functions to be 


/».(*) = 


• (°) 
2 


T.<°> , 

2 


We define the operators A it using (3.19) as follows: 


Au — r u ( 0 ) + T ' u ( 0 ) * 


1 1 


_L.Vff ZgL'dtdt. 

4*11 r <j 


U 00 


(3.20) 
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Finally we take T to be the property that the following 
inequalities are satisfied: 

r ks (t)>fUt). (3.2D 

The conditions of Theorem 3.2 are satisfied since, from the 
assumption that (3.21) are satisfied in the interval 0 ^ t < r, 
it follows that 


i i 


1 

a k$ 


o 00 V 


t t 


<>0 0 


Therefore, by Theorem 3.2, T is an inductive property 
for those sytems of gravitating bodies for which the 
inequalities 


2 ~ 1_ 2 1 



numj 


dt dt > 0. (3.22) 


are satisfied for all f^O. These inequalities should be 
considered as defining the appropriate values of r,^(0) 
and r'ijiO). 

Evaluating the integral on the left-hand side of (3.22) 
we find 
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t t 

mm 4 

- — r± dt dt 

r <j(°) , r o(°M 

2^2/ 

4/n»my ^ km^mj 

~ r (J (0) 4 (0) ,J(0) 

Substituting the value of the integral in (3.22), we obtain 
after simplification 


In 1 + 


'o-(O) 

'</«» 



1 "V»_4nnwy ) 

1 2 ^ 1 ? ^(0)^(0)j 



The second and third terms on the left-hand side are 
positive for every value of t>0, provided that t ij (0)>0. 
Thus, (3.23) will be satisfied for all positive t if the co¬ 
efficient of t in the first term on the left-hand side is 
positive, i.e., if 

2 r ij (°> 

When these conditions are satisfied, inequality (3.21) 
will hold for all positive t, provided (3.21) holds for t=0. 
But then (3.21) is satisfied at /=0 since 
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/,.<oj=^a. 

We can thus formulate the following theorem. 

Theorem 3.4. If the conditions 

r' 0 ( 0 )> 0 , 

r 'ij (°) _i_ ^ . q (3.24) 

are satisfied at the initial time t = 0 for a system of gravitating 
bodies , then the distances between the bodies tend to 
infinity as t-+ . 

5. The criterion (3.24) just obtained is rather simple, 
but, nevertheless, slightly cumbersome. However, by 
weakening it, we can reduce the criterion to a simpler 
and more transparent form. 

Let us introduce the notation 

p(/) = mjn [r {J (t)}, 

-. 

min ynj 

It is not difficult to see that conditions (3.24) hold if 

a(0)>0, 

s>(0)>Hl, (3.25) 

w 'p(U) 

Thus, if conditions (3.25) are satisfied, the distances r i5 
between the bodies tend to infinity as . 



METHOD OF CONTINUOUS INDUCTION 


41 


We might note that conditions (3.25) have exactly the 
same analytic form as those which were obtained in 
Chapter 2 using dimensional analysis. 

6. Let us now consider a system of gravitating bodies 
Pot Pi» . .., P n — i, the motion of which we shall describe 
using Jacobi coordinates. Set p = p«— 1 » where —i de¬ 
notes the distance of P n ~.i from the center of mass of the 
remaining bodies. Choose a positive number a^l, and 
consider the sphere 5(a) with center at the point P n -1 
and with radius P = ap. Thus, the center of the sphere 
moves together with the point P„_i, and its radius varies 
with time. We shall say that an a -approach occurs be¬ 
tween the body P n -1 and some other body P,- at time t 
if at this time Pi lies inside or on the boundary of the 
sphere 5(a). 

Differentiating the identity 


P 2 = 55—1 + 1 + 1 


twice with respect to the time, and carrying out the same 
kind of manipulations as was done for foi in Section 4, 
we obtain 






dU 
ol » 


(3.26) 


where l is the direction of the line joining the center of 
mass of the bodies Po, Pi> • • ■» P 2 to th e point Pn—i* 
Since in differentiation with respect to the direction l 
only the distances between P n -i and the remaining 
bodies vary, we have that 


H-I 


dJJ_ 

fit 





r 


2 


w-1 


1 

dl 
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From simple geometrical considerations, it follows that 



If it is assumed that no a-approach between P n _i and 
the remaining bodies of the system occurs when 0, then 

r i>n ^ { >ap. 

Using this to strengthen the previous inequality, we obtain 

(3.28) 

where 

Finally, if we integrate (3.28) twice with respect to time 
from / = 0 to />0, we obtain the integral inequality 

P<0>P(0) + p'(0)*-£ (f * (3.29) 

o o P 
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Theorem 3.5. If the initial conditions 

P '(0)>0, 

f''(0)-Hi)> 0 ' (3 ' 30) 

hold att = Q for a system of gravitating bodies P 0 , Pi,..., P„— 1 
and if no a-approach of P n —i to the other bodies occurs for 
any value of- f ^ 0, then the distance of P n —\ from each of the 
other bodies in the system becomes infinite as f-> «. 

Proof: The integral inequality (3.29) is of fhe same 
type as those in the system (3.19). Applying Theorem 
3.2 to it and repeating the argument that was used in the 
proof of Theorem 3.4, we obtain a proof of the stated 
theorem. 

7. We now consider some applications of the second 
Theorem on the Inductiveness of Inequalities. As a first 
example, it will be constructive to consider one of the 
cases already treated with the help of the first Theorem 
on the Inductiveness of Inequalities. Let us turn again 
to Theorem 3.5, and let us prove it in a slightly stronger 
form. In so doing, we obtain some results which we shall 
find useful below. 

Theorem 3.5a. If the initial conditions 
p'( 0)>0, 


hold at t =0 for a system of gravitating bodies P o, Pi, P *—u 
and if no a-approach of P*_i to the other bodies occurs for 
any value of t^ 0, then the ineQuclities 
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(3.31) 

p(l)>f(0) + t^l, (3.32) 

hold for all positive t, and, therefore, p (/)—>■ <» as t—r oo. 
This implies that the distance of P K _! from each of the other 
bodies in the system becomes infinite as t-y co . 

To prove the theorem, we start with inequality (3.28). 
Integrating it with respect to time from t = 0 to t>0, we 
obtain 

o Y 

Applying the parallel version of Theorem 3.3, we set 

*W=p(0. 

and we define the operator A by means of the equation 

II 

For the comparison functions, we take 

tw=qa. 

/(0=p(0)+E^!(. 


It is first necessary to show that the conditions of 
Theorem 3.3 are satisfied for the case in question. With 
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our choice of 8(t), w(/), <p{t)> and /(/), it is obvious that 
(3.9) holds. The operator A also satisfies the required 
conditions since from the asumption that p(t)>f(t) in 
the interval 0^<r, it follows that the inequality 




t 



is satisfied for these same values of t. Therefore, by 
Theorem 3.3, inequality (3.31), as well as inequality (3.32) 
which follows from it, are inductive properties for systems 
of gravitating bodies for which the inequality 


¥■ 



> 0 , 


holds for all values of 0. After integration, this last 
inequality becomes 


Am i_2A/_. 

2 ' *'<'» + 4 (0)+ £ip,) > • 

It is easy to see that when the conditions of the theorem 
are satisfied, this inequality holds for all values of t ^O. 
Therefore, under the conditions of the theorem, (3.31) 
and (3.32) are inductive properties, and since they hold 
at the initial time /=0, they are also valid for all 0. 
The theorem is proved. 

8. Let us now consider some further applications of 
the second Theorem on the Inductiveness of Inequalities. 

Suppose we have a system of gravitating bodies 
•Po, Pi, . . P n —i with masses wo, Wi, . .To de- 
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scribe their motion, we make use of equations (1.12). 
Integrating these equations with respect to time from t=0 
to />0, we obtain 



From this, it follows that 



With a view towards applying Theorem 3.3, we set 

“.>(') = '■’,( 0 . 


and we define the operators A ({a,-*}),{«}) by means of 
the equations 




(3.34) 
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We can thus consider (3.33) to be operator inequalities 
satisfied by the system of gravitating bodies. 

We now assume the initial conditions are such that 


**<0)“2r*(0»0. 

*4(0)-r;*(0)>0, 


(3.35) 


and we define the comparison functions by means of the 
equations 


?«w=4|><0)+<*<o)]. 

/*(,) = [ |v tt (0)-|r; t (0)]V+ 

+ 2r tt (0) [| (0) -1 (0)] t + H* (0). 

Finally, T will be considered to be the property that the 
inequalities 

K,<0- r ; 1 <°)|<T[r 1 <0)+'',(0>]. 

hold, where v ik ( 0)»|r'^(0)|. 

To prove the inductiveness of I\ we require the 
following lemma (G.A. Merman, 1953a). 

Lemma 3.1. If the vector r (t) satisfies the inequality 

| r'— r | (0)<€ 

in the interval 0^/<r, then it satisfies the inequality 
r 2 > [t>(0)— e] 2 t 2 + 2r(0) [^(0)—e] Z+r^O), 
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for these same values of /, where v(Q) = |r'(0)|. 

Proof: From the integral identity 

t 

r~ r{0) + r' (0) / -j- J [r' — r' (0)] dt 
we obtain 

r > I r (0) + r' (0) /1 -1 J (r' - r' (0)J dt | , 

o 

Hence, by the hypothesis of the lemma we have 
Olr(0) + r'(0)*|-e* f 

or 


r2 > l r (0) + r ; (0) t] 2 — 2e | r (0) r' (0) /1 ^ —|— eV, 
which proves the lemma when the inequality 
|r(0) + r'(0)/|<r(0) + i;(0)/. 


is applied. 

Using this lemma, we now show that if the property 
T, i.e., the inequalities 

| 4 - 4(0) KjK (0) -f 4 (0)], (3.36) 

hold in the interval 0^i<r, then the inequalities 

4 > i 4 »•« C») — j 4 (0) f -f 

^ 2r ‘» (°) [ (°) - J ”«(0)) t + r% (0). (3.37) 

will also hold for these same values of t. 
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To prove this fact, we apply Lemma 3.1 for which it 
suffices to set 




and to carry out some simple manipulations. 

From the expression (3.34) which defines the operators 
Aij } it immediately follows that if it is assumed that 
0 * n the interval 0^/<r, then the inequalities 



will hold for these same values of t. 

The results obtained show that the conditions of the 
second Theorem on the Inductiveness of Inequalities are 
satisfied by the case under consideration. Therefore, the 
inequalities (3.36), as well as the inequalities (3.37) which 
follow from them, are inductive properties for systems of 
gravitating bodies which satisfy the inequalities 


n—1 


-f£+2«4i&+ 




n— i r 

4- 2M JZ — 7 ^rt(°)4-r.*(0)] < 0. 

j=\ n * 


(3.39) 
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for all values of / ^ 0 

Moreover, (3.36) and (3.37) are automatically sat¬ 
isfied at the initial time / = 0; therefore, a system of grav¬ 
itating bodies satisfying condition (3.39) will satisfy (3.36) 
and (3.37) for all values of 0. Thus, the distances be¬ 
tween the bodies tend to infinity as t—> °°. 

It remains for us to investigate conditions (3.39). 

All the integrals which appear in (3.39) are of the 
same type* namely, for any values of i and j we have 

t t 

| _f_ dt _ 

lij J «ijt* *+ bijt -f Cij ’ (3.40) 

u o 


where 


= [ f V H (°) — T r ‘J (°) J. (3.41) 

K = [ T r '<J (°) — J V *J ( ( >)] 2r (J (0), (3.42) 

r U = r i ij(0). (3.43) 

It is well known that the properties and the evaluation 
of any integral of the form (3.40) essentially depend on 
the sign of the discriminant 

Using the values of a ijt b ijt and c, j given by (3.41) (3.43), 
as well as conditions (3.3S), we find after some simple 
computations that 
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\j =- 4 ^ <°> K; (°) “ r J (°)J > °- O’- 44 ) 

For A,v>0and a,;>0, the quadratic form a, ; / 2 +•/;, ,/+- 
Cj-y is ix)sitive for all values of t. Therefore, 

t CO 

I aijt* H- bijt + C{/ ^ 1 4- 5iyl + c*/ (3.45) 

0 0 

Furthermore, when A f J >0 and, as in our case, a l; >0, we 
have 


CO 


f 


dt 

+ bijt -f Cij 



arc tan 


4*4. 


Since 0, as follows from equation (3.42) and condi¬ 
tions (3.35), from the last equation we obtain the estimate 

on 

r lit ^ IT. 

J b,jt + Cij ^ v to, 7 ' (3.46) 

II 1 


Finally, from (3.40), (3.44), (3.45), and (3.46), we find 
that the inequalites 



*2 7T 

<«) K»?, <‘>) - K-J (°) 


(3.47) 


0 
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hold for any values of i and j and for all 0. From this 
last result, we immediately obtain the estimate 


t n—1 * n^\ f 


J= 0 




<T 

V3 


n(.U — »),) (.U — J-M1-1 
(0) )/ (i ) -r- k {0) 


(3.4S) 


where M- YA~h m i’ 

On the basis of (3.48), we can state that if the 
inequalities 


r. (Af — nu) (M — mj () _ 

* ^K 

-T['.< 0 )+ r «(°:)]<“• 

hold at the initial time / = 0, then (3.39) and all other 
consequences of it will hold for all values of t^O. 

To conclude our investigation, we can formulate the 
following theorem. 

Theorem 3.6. If the initial conditions 


(0) 2r' k (0) > U, 


'^(0)-4(0)>0, 


‘■*(0)+'•*(“)]K^(0) -r;*(0)> 3v3 " (it ' 
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are satisfied at t = 0 by a system of gravitating bodies, then 
the distances r ik between the bodies tend to infinity as t-+a . 

We shall now make some comments apropos of the 
theorem just proved. 

The criterion we have obtained for the complete 
scattering of a system of bodies as <--><» requires each of 
the velocities v ik (0) to be sufficiently great in comparison 
with both the reciprocal of the distance r ik ( 0) and the 
radial velocity r',*(0). The first requirement is intrinsic 
to the problem, whereas the second is a consequence of 
the mathematical reasoning used in the derivation of 
the criterion. 

The theorem imposes the fewest conditions on the 
initial state of the system when r',*(0) = 0. We note, in 
passing, that it is just in this very case that the criterion of 
Theorem 3.4 fails to answer the basic question. When 

(0) = 0, the first and second conditions of Theorem 3.6 
are automatically satisfied and the third simply becomes 


_* tO\>i 3rtr«(Af-.!*)(*-**) 

** (0) > ..r«*(0)- 

If the relative and radial velocities are equal, namely, 
if v ik (0) = r' ik (0), then Theorem 3.6 fails to produce an 
answer. However, in this case, Theorem 3.4 is appli¬ 
cable and the criterion established there is the most suit¬ 
able one. Thus Theorems 3.4 and 3.6 complement one 
another. 

Finally, we note that the criterion obtained above may 
be regarded as a weakened variant of the criterion for 
hyperbolic motion formulated by G.A. Merman (see the 
reference 1953a) for the particular case of the three-body 
problem. We would like to remark on this variant in the 
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following connection. The criterion of Merman is not 
very transparent, and its form does not permit us to tell 
whether the set of initial conditions satisfying this criterion 
is not empty. In our case, this is obviously so. More¬ 
over, it is not clear under what conditions Merman’s 
criterion is applicable; however, they cannot differ essen¬ 
tially from those under which Theorem 3.6 holds. 

9. Let us now consider a system of gravitating bodies 
.Po, Pi, . . ., P n — i- We describe the motion of the system 
using Jacobi coordinates. Set p = p,_,, where p„_i de¬ 
notes the distance of P„_i from the center of mass of the 
remaining bodies. Let p denote the vector joining the 

center of mass of the system P 0 , P,.P„_ 2 to the body 

P„_i. Then 


P—*’«-i H” j 7 i«-i + (3.49) 

where i, j, and k are unit vectors directed along the 
£»— i, Vn—i, and f n _i axes, respectively. Differentiating 
(3.49) twice with respect to time, we obtain 


l^iC-i+KU + kw 


or, by the equations of motion (1.7), 


P*= 




(3.50) 


Now, using the relations 
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»-* 


dU _ 

— i 

^ r *»-j 

<Jr<, w -i 


N—2 

y^-i 
<=0 i 

w—l 

l 

T 

*v.i 

_d£__ 

»-2 

r *"-i 

<?ri. n— i 




we can write (3.50) in the form 



Integrating this equation with respect to time from 2 = 0 
to 2>0, we obtain 

p'—P'(0)= 

1 f | . dr <, »-l i —l 1 „ 


_ _1_ f "V 1 J <■ «-l I : »-l i t z± "- 1 

d5 -i d \-i ,K »- 1 


Using the fact that 


I dr t.n-l I . dr <, n-l | . dr (, .- 1 

*_ *W_! 


we find from this that 


ip'-rwK si— J2 : ^ ===i ‘"- 

0 4=0 *—1 


(3.51) 
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If it is assumed that no a-approach of JP»— 1 to the 
other bodies occurs for t^O (the definition of a-approach 
was given in Section 4 of this chapter), then 

1 > ap, 

Using this to strengthen inequality (3.49), we then obtain 

Ip'—|f<0)l<^rj jt. <*•“> 

0 


where 


M* — 




—5 




The integral inequality (3.52) is of the same type as those 
in the system of integral inequalities (3.33) considered in 
the preceding section. Applying the second Theorem on 
the Inductiveness of Inequalities to (3.52) and repeating 
the argument used to prove Theorem 3.6, we arrive at the 
following theorem. 

Theorem 3.7. If the initial conditions 

»(0)-2p'(0)>0. 


i' 2 (0) — p' 2 (0) > 0, 


l '(0)+f/ (0)1 Vr*(0)-:/*(0); 


:i v’3 r. 
( 0 ) 


are satisfied at t=Q by a system of gravitating bodies , where 

•«=vc,+c,+c.. 
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and if no a -approach of P n —i to the other bodies of the system 
occurs for any value of t^ 0, then p (/)-> od as t-+ co. There¬ 
fore , /fce distance of P n _! /row eacA 0 / /fte other bodies will 
become infinite as /-> co . 

10. We shall say that a system of gravitating particles 
is completely dissipative as 00 if all the distances be¬ 
tween the particles increase indefinitely as 00 . 

Theorem 3.8. Jf aZ the initial moment of time /= 0, the 
conditions 


P'(0)>0, 

<>'•«•)> < 3 - 53 > 

ore satisfied by a system of gravitating bodies 1 P 0l Pi,..., P*_i, 
wAore if is //re totaZ energy of the system , and if //tore occurs 
no a-approach of Pn—i to the other bodies for any value of 
t^ 0, Zften ZAe distance of P n —\ from each of the other bodies 
increases indefinitely as oo, but the system of bodies Po, 
Pi, . . Pnr-i on the whole is not completely dissipative . 

That the distance of P*_i from each of the other bodies 
in the system increases indefinitely as q° follows im¬ 
mediately from Theorem 3.5a. 

We must now show that the system of bodies in quest- 
tion cannot be completely dissipative as /->-». 

From the energy integral 

.=-1 


we have 
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V > 1 (Cl + C. + C.) 


and using the obvious inequality 

«*_i + Ci +Ci > Cco. 


we find that 


V)-H. 


From Theorem 3.5a it follows that p'(t)> J p'(0) for all 
values of / ^ 0. Therefore, 

Hence by (3.53) we have 

U> 0 

and, therefore, U does not approach zero. This means 
the system of gravitating bodies P 0 , Pi, . . P„_i cannot 
be completely dissipative as t-^oo. 

11- In proving some of the theorems, namely,Theorems 
3.5, 3.5a, 3.7, and 3.8, we assumed there occurred no 
a-approach of P w __i (whose motion we were studying) to 
the other bodies Po, Pi, . . ., P«— 2 * It seems to us that 
the presence of such a condition in the problem is due to 
its intrinsic nature. 

In fact, consider the aggregate of the gravitating bodies 
Po, Pj,..., P n _i, and suppose that P w _i is located inside 
the aggregate. If the initial velocity of P w _ x is sufficiently 
great, and if while moving inside the aggregate P n _i does 
not approach the other bodies closely, then P n ^ x will leave 
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the aggregate and go off to infinity. In other words, its 
distance from each of the other bodies will become in¬ 
finite as /--> oo. 

However, if P n —\ comes close to the other bodies, this 
could result in such an exchange of energy between the 
bodies in the system, that P w — 1 , by having its velocity 
changed, could not leave the aggregate in spite of the 
fact its initial velocity was sufficient to have this happen 
under the other conditions. 

However, it is necessary to observe still one other 
aspect of the problem. It could happen that two of the 
bodies P { and P^i are ejected from the aggregate, and in 
moving an infinite distance away from the other bodies 
of the system, are within an a-proximity of one another 
or else arrive at this state after a certain time has elapsed. 
In this case, the theorems we have proved fail to solve 
our problem. 

Thus, the criteria formulated in Theorems 3.5, 3.5a, 
3.7, and 3.8 yield a solution to the problem only if the 
body in question experiences no subsequent a-approaches 
to the other bodies. Since these theorems do not indicate 
under what conditions any a-approach is excluded, they 
merely have qualitative significance. Nevertheless, they 
can be used quantitatively in some problems in the dy¬ 
namics of stellar systems. For example, by knowing the 
parameters of a stellar system, it is possible to estimate the 
probability of an a-approach of a given body to any of 
the other bodies and, therefore, to estimate the prob¬ 
ability of occurrence of the type of motions we have been 
considering, under given conditions on the initial state 
of the system. 

Finally, we note that the condition that there be no 
a-approaches can be omitted in certain instances since it 
is automatically satisfied. Namely, this is so in the 
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application of Theorem 3.8 to the particular case of the 
three-body problem. It was in this way that the author 
first obtained criteria for hyperbolic-elliptic motions in 
the three-body problem. Later on, G. A. Merman, using 
our method of continuous induction, gave several inter¬ 
esting theorems containing more exact criteria for hyper¬ 
bolic-elliptic motions in the three-body problem. 



CHAPTER 4 

Method of Invariant Measure 

1. The application of the method of invariant measure 
to mechanical problems and, in particular, to the man y, 
body problem is based on the phase interpretation of 
mechanical systems. 

The phase interpretation of a system of « gravitating 
bodies is as follows. In Chapter 1, we saw that the 
second-order differential equations of motion can be re¬ 
placed by an appropriate number of first-order equations, 
the general form of which may be written as 

• • • i *,), 

. (4.1) 


The state of the system of bodies is then completely 
defined at each moment of time t by the values of v real 
continuous functions of /, *i, x 2} . . . , x ¥ . We call each 
possible set of values of x h x 2t . . . , x ¥ a phase of the 
dynamical system. While regarding the collection of 
numbers x %, as <the phase of the dynamical 

system, we shall, at the same time, think of it as a point in 
p-dimensional euclidean space which we refer to as 
phase space . 

We now choose some moment of time as our initial 
time which we designate to be ^ = 0. Each point p£E\ 
with the exception of those which are singular points of 


61 




62 QUALITATIVE METHODS IN THE MANY BODY PROBLEM 


the given system, is a possible phase of the dynamical 
system at time t = 0. The phase of the system at any time 
t other than the initial time will in general be different 
from the initial phase. We let denote the point of 
Er that represents Jthis new phase, thus f(p, 0) = p. As the 
time varies positively, will be defined in some in¬ 
stances for all positive values of t , whereas in other 
instances, there will be a r>0 such that fip>t) will be 
defined only for all positive t<r. A similar situation 
will prevail when the time varies negatively. 

As the time varies, the phase Xi, x ¥ will change 

in accordance with the laws inherent in the motion of the 
mechanical system. Corresponding to these changes of 
phase, the point/(/>,/) will move in a given way in E v and 
will describe a certain trajectory. 

By simultaneously considering all the points of phase 
space and their motions, we obtain the picture of the flow 
of the “phase-point fluid.” This picture describes the 
aggregate of possible motions of the dynamical system. 
By applying the method of invariant measure, we shall 
be able to characterize the motion of the particles in 
physical space in terms of the motion of the points of 
multidimensional phase space. 

2. For each pair of points of the space E? 

P ~( x ii 

P" = ( X U X t, 

we introduce the nonnegative number 

p(p',p")=V2K-<)*’ 
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called the distance between p ' and p". The general 
properties of the distance in E“ are exactly the same as 
those in three-dimensional space: 

(1) (Identity property) p{p',p") = 0 if and only if 

P • -P . 

(2) (Symmetry property) p(/>',£") = p (p",p '); 

(3) (Triangle property) For any three points p\p" >p ,n 
of E", we have 

p(P'> P")+P(P", p m )^?{p\p m \ 

Let p be a point and A an arbitrary set of E\ The 
nonnegative number 

'4P, A) = inf p(p, q) 

q£A 

is called the distance between the point p and the set A. 

Using the notion of distance, it is easy to introduce into 
E v the notions of the neighborhood of a point and the 
neighborhood of a set. If p is any point of E v and e is 
an arbitrary positive number, then the set of all points p ' 
for which p(p,p')<e is called a spherical neighborhood 
of p of radius e and is denoted by S(p , e). If A is any set 
of E y , then the set 

S(A,t)= U S(p, e) 

is called a spherical neighborhood of A of radius €. 

Although it is impossible to give a graphical representa¬ 
tion of the motion of a point and its trajectory, or the 
motion of a set, in phase space, nevertheless, both the 
above interpretation with its geometrical terminology and 
the analogy of the properties of the space E v to those of 
ordinary three-dimensional space considerably facilitate 
the discussion. 
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Since distances and neighborhoods are characterized in 
y-dimensional space by properties analogous to those of 
distances and neighborhoods in ordinary spaces, the theory 
of point sets and the theory of measure can be suitably 
generalized to the space E v . A presentation of these 
theories and, in particular, the basic concepts that we 
shall use may be found in books on the theory of sets and 
the theory of functions of a real variable. 

3. The general properties of the phase motions j(p,t) 
will, of course, depend on the character of the functions 
X{(x i,... , x v ) that appear on the right-hand side of (4.1). 
For a wide class of equations of the form (4.1), to which 
the equations of motion of a system of gravitating bodies 
also belong, it can be shown that the phase motions have 
the following properties: 

(1) Uniqueness property —every point p of phase space 
(with the exception of singular points) defines a unique 
motion of f(p,t) in an interval of time of one of the follow¬ 
ing four types: 

— 00 </<+ °° , 

-OO <t<t'\ 

t'<t<+ CD, 

Thus, through each point of phase space passes only 
one phase trajectory. 

(2) Stationarity property —in the interval in which 
f{p,t) is defined, the condition 

f(f(Pi ^i)» h)=J(pi 

is satisfied. 

(3) Continuous dependence of the motion on time —for 
every / in the interval of time in which f(p,t) is defined and 
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for any € > 0 t we can find a 5 > 0 such that 

P(f(p, 0 . /(/>/))<€ 


for all /' such that 


I t-t' | <5; 

(4) Continuous dependence of the motion on the initial 
point —if the motion of f(p t t) is defined over the interval 
[tiA]> then for e >0 we can choose a positive 8 sufficiently 
small so that for any point q for which f(q,t) is defined and 
satisfies the condition 

POXff. h ))<«, 

the motion of f(p t t) exists in the interval [h,h] and satisfies 
there the condition 

P(/(?> <))<«• 

4. Of great importance below will be the concept of the 
measure of a set of points of phase space. We shall use the 
measure normally defined in multidimensional euclidean 
space by the direct generalization of the concept of 
measure in three-dimensional space. This measure is 
only meaningful for bounded sets. However, along with 
bounded sets, we are going to consider unbounded meas- 
surable sets both of finite and of infinite measure. To 
our end, it suffices to define the measure of an unbounded 
set as follows. Let 22i, J3 k, . . . , jR»» • • - be an unbounded 
monotonic increasing sequence of positive numbers. Let 
S n denote the sphere of radius R* with center at the origin. 
We call an unbounded set A measurable if each of the sets 
A OSn is measurable. The sequence meas(-<4 OS’*), n = 1,2, 
. . . , either tends to a limit or becomes unbounded as 
w-»*oo. If the limit exists, we then define 
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meas A -lim meas (^HS*). 

and in the contrary case, we set 
mes A = oo 

Clearly, when this definition is applied to a bounded set, 
we obtain its usual measure. 

The great importance of the concept of measure is in 
the fact that Liouville's Theorem on Invariant Measure 
is applicable to the motion of a system of gravitating 
bodies in phase space. The invariance of measure in phase 
space is the basis that allows us to ascertain some of the 
very important properties of dynamical systems. 

Theorem 4.1. (Liouville’s Theorem) If the equations 

§ = *(*1 .*.), 


d*, 

dt 


—X,(x u 




(4.1) 


are such that 


V 



(4.2) 


and if A is a measurable set of E* of finite measure such that 
for any point p € A the motion of f(p t t) is defined for all 
t^O(t^O), then 


meas f(A, /) = meas^4 


for all ^0(/^0), i.e., the measure is invariant vihen 
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Proof: Let A be a measurable set of E" of finite measure, 
where for any p€A, the motion of f(p,t) is defined for all 


mea s f (A, t)= J dxi ... dx,. (4.3) 

*A,i) 

Consider the solution of (4.1) 

fll, . • i | 

. (4.4) 

~ x* ( t ; fli, •••* flv), 
satisfying the initial condition 


= a* ; 


for / = U. 


where (a lf . . . , a v ) is any point of the set A. 

Let us introduce new variables of integration xi t ..., 
defined by equations analogous to (4.4): 


X\ — X\ i X \» ■ * • * x v )» 


Xv —— Xv (/j X\ , • • • | Xv), 


(4.5) 


in which / is regarded as a parameter. In terms of the new 
variables, (4.3) becomes 


meas f(A, t)= ... dx,, (4.6) 

4 


where D is the Jacobian 
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Q — d <*i» •••■ *»> 

~i*i .<r • 

Differentiating (4.6) with respect to the parameter t, we 
obtain 

. . di„ 

A 

Then going back to the old variables x u ..., X<1 we find that 

iwas/M, 0= J -±-£dx, ... dx.. (4.7) 

By the rule for differentiating a determinant, 


where 

0 < = 

or by (4.1) 



do ^ , 

It “ > L) >- 



r. ^4% 

, *, +1 . 

1 ' '» 

d(x h 

.... J.-J, 27 , 2-,. + lt . 

* * * ) 

*(*„ . 

*<—l, Z4l, .. 

1 • | 

<>(* 1 , . 

• •. i ( , i <+]l ,. 



(4.8) 


However, since 
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££i j£u. 

<*** *(*», 


By then noting that 


J <-n *>» A-t-it *“t *J 

*<» *<+i* *”i 


Ifgii • • •» y <-i«*»• -ti-fi. •.x„) __ f for ik«) 

'MA* ■ • •'A-t. Ai A,it....#,) j 0 {of 


we find 


"< 



Substituting this expression for D, in (4.8), we obtain 
— — /> \i£i 

<#< ~ 2L ax, * 

<«i 


and hence 


1 <)/> _<*A ( 

TT at ~ * 

<=i 

Using this last result, we can express formula (4.7) in the 
following form: 

/• * 

^jmeas/M, t) = l 24iT <**!••• (4-9) 

J <*■! 

/(ii« 

This is the remarkable formula of Liouville. 

We now proceed to prove Liouville’s Theorem. From 
formula (4.9) and condition (4.2), we have 
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— mens /(A,t) = 0, 

Integrating this with respect t from J = 0 up to />0 and 
taking into account that f(A } 0)=A, we find that 

meas f (A, /) = meas A. 

The theorem is proved. 

It is necessary to emphasize the importance of the 
requirement that for any point p€A ,the motion of f(p,t) is 
defined for all i^0(/^0). If this requirement is not 
satisfied, the set f(A,t) and its measure may not be 
meaningful for all values of 0(/^ 0). This can even be 
true for differential equations of a very simple form as the 
following examples show. 

Example 1. Consider the motion defined by the 
equations 


dx 

~d7 



(4.10) 


Phase space for this system is the euclidean plane E 2 . 
From (4.10), we find that 



y = Ce-\ 


where C is the constant of integration. Thus the trajec¬ 
tories of the phase motions are exponential curves. Let 
us choose an arbitrary point p in the plane with coordinates 
xo t y 0 , and let us consider the motion of f(p,t). Integrating 
the first of equations (4.10) subject to die initial condition: 
x=x 0 when / = 0, we find that 

x = —In (e~~ x o — t). 
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From this result it is easy to see that as *o. 

Therefore, the phase point f(p,t) goes off to infinity in a 
finite interval of time, and for all o, the motion of 
f(p,t) is not defined. This is true for any choice of the 
initial point p€E 2 . 

Example 2. Consider the motion defined by the 
equations 




< { y 

dt 


— — \'2~ xye 2 


(4.11) 


Phase space for this system is again the euclidean plane E 2 . 
From (4.11), we find that 


V = £ e 


where C is the constant of integration. Thus the trajec¬ 
tories of the phase motions are normal curves. Let us 
choose an arbitrary point p in the plane with coordinates 
xo,yo , and let us consider the motion of Integrating 

the first of equations (4.11) subject to the initial condition: 
x~xq when f=0, we find that 



*0 


From this we see that 
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X 4- oo as 



—oo a ? 3 

*— —« as t - -L- f r" rf:r>—1. 

Tlius, the phase point goes off to infinity in a finite 
interval of time when the time varies both positively and 
negatively. In this example, every motion is defined in 
an interval of the form 


where 0—1 and <"^1. 

5. The main object of study in general dynamics are 
the invariant seta of phase space. A set A is called positive 
(i negative ) invariant or invariant as t-> ® (/->— oo ) if for 

every point p£A, f(p,l) i s defined for all /^O(f^O), and 
for these times 


fiPMA. 

A set 4 is called simply invariant if it is invariant both 
as t-+c o and as t->— « . It is clear that the union or in¬ 
tersection of any collection of sets that are invariant either 
39 . °r as oo , as well as the difference of two in¬ 
variant sets, is also invariant as or as t-+— oo. 

«?***•** of aI1 P° ints PZE’ for which 
Ae motion of MO is defined for all t> 0 (l$0), and let 
w U4 1 denote this set. Obviously, S + (S“) is a 
maximal set, invariant as /->» (<-*_ „ ). we rail it th* 
regular kernel of the dynamical system for tU- oo (*->_ ») 
We call the set } ‘ 
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s=s + ns- (4.i2) 

simply the regular kernel of the dynamical system. 

Theorem 4.2. E + , E“, and E are sets of type G Sy 
i.e. } sets which are the intersection of a countable number of 
open sets. 

We first prove the theorem for the setE + . Let 
Ti, T 2 , . . ., t w , ... 

be an unbounded monotonic increasing sequence of pos¬ 
itive numbers as n -+Let E t denote the set of all 
points p of phase space for which the motion of f(p t t) is 
defined in the interval [0 t r h ]. If pZEt, then due to the 
continuous dependence of the phase motion on the initial 
point, we may choose a positive number 8 k such that for 
every point q , with g€S(£,5*), the motion of f(q,t) is 
defined for all 2£[0 ,t*]. However, this means that 
S(q, , and therefore Ej£ is an open set. Obviously, 

CZ oo 

X-l 


from which it follows that E + is a set of type G s and is 
therefore measurable. 

The theorem is proved in an analogous way for S 

The validity of the theorem for the set E follows 
immediately from (4.12). 

6 . Let us take some set A of E v at the initial moment 
of time. Each point p€A will move in its own way in E ", 
and at some other moment of time, all the points p of A 
will occupy new positions f(p,t) and form a new set; we 
denote this set by f (A } t). 

Let p and A be such that f(p,t) a.ndf(A f t) are defined 
for all 0(t ^0). From the definition of/(/>,/) and/(-4,/). 
it follows" that if p€ A, then f(p,t)€f(A,t) for every 
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/^0(/^0). Conversely, if f(p,t)£f(A,t) for some 
^0(^0), then p€A. 

From this remark, we immediately have the following 
simple facts. If Ai and A 2 are sets such that Aif^Az^O, 
and if f(Aut) and f(A^t) are defined for all /^0(/^0), 
then for these values of t we have 

/(i4i,0rWt,0«0. 

On the other hand, if Ai and A 2 are sets such that A iDA* — 
0, and if f(Ait) and f(A 2} t) are defined for all /^0(/^0) 
then for these values of t we have 

f(Au t)r\f(A 2} t)y£ 0 . 

A set A is called recurrent as /-> 00 (/-->— 00 ) if for any 
7>0, we can find a /'> T such that 

Anf(Aj')^ 0 (AnKA-n* o). 

Conversely, the set A is not recurrent as /-><» (/ ■>— 00 ) 
if we can find a r > 0 such that 

AClf(A % t)-a (An/(A,-t)-0) 

for all t>T. 

Theorem 4.3. If a set A is recurrent as /->» (/->— 00 ), 
and if /(/>,/) is defined for all /<0(/>0) and for any point 
p£ A, then A is recurrent as /->— cd (/-> » ). 

If A is not recurrent as /-> co(->— 00 ), and # /(/>,/) 
defined for all /< 0 (/> 0 ) and for any point pO.A , /Aen 
^4 is no/ recurrent as /->— ® (/-> 00). 

Proof: Let the set ^4 be recurrent as /~>oo. Then, for 
any positive T no matter how large, we can find a /'> T 
such that 


Anf(A,t)*0. 
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We then have 


f(A,—t)nA?£ 0 , 

and this implies that A is recurrent as t ~>— ® . 

For a set A which is recurrent as /->— ®, the theorem 
is proved analogously. 

The second part of the theorem is easily shown by 
arguing by contradiction. 

Theorem 4.4. If a set A is such that f(A,t ) is defined 
for all real values of t, and if we can find a number r>0 
such that ADf(A,t) = 0 or A(~)f (A—1)=0 for all t^T, 
then the sets 

.... f(A-nr),... ,f(A ,-t),A ,f(A ,r),... ,f(A , nr ),... 
are pairwise disjoint. 

We prove the theorem by contradiction, and we assume 
that we can find two integers i and j, with i<j, such that 

f(A,iT)C\f(A,jT)^0. (4.13) 

Consider the sets A and f{A,{j-i)r)\ since j-i ^ 1, we have 
by assumption 


ADf(A,(j-i)T) = (). 


and, therefore, 

f(A, ir)C\f (f{A ,(j — i)r ), t'r) = 0, 


However, 

f(f(A , (j — i)r ), ir)=f(A, jr). 


Thus 


f(A,ir)nf(A,jT) = 0, 
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and this contradicts (4.13), and so the theorem is proved. 

The theorem is proved analogously for the case where 
Ap\f(A ,—0 = 0 for all t^r. 

7. A point p is called recurrent as t-+ oo (/-->— oo ) if for 
any positive T no matter how large and any € >0, there 
exists a time t f >T(t'< — T) such that 

f(P,0£S(p . e). 

A point p is called simply recurrent if it is recurrent both 
when t-+ oo and when —« . 

Theorem 4.5. If the point p is recurrent as t—>° ° 

0 t — co), then for any ti>0(ti<0), the point f(p } t i) is 
also recurrent as /-->• oo (/->— oo ). 

We shall only prove the theorem for a point p which is 
recurrent as /->co . The case where the point is recurrent 
as /->—oo is treated analogously. 

Let us take arbitrary positive numbers € and 7\ and 
let us choose 5 = 6 (e ,T) sufficiently small so that 

f(S(P< 5)/0CS(/(A *i),€). 

This is possible because of the continuous dependence of 
the motion on the initial point. Since the point p is 
recurrent as oo , we can find a t'>T such that 

f(P/)£S(P,8). 

If Q is a point of the set A y then by the definition of 
and/C4,/) it follows that f(q y t)£f(A,t). Using this fact, 
we can write 

MP> a h) € f(S(p, 5),t 1 )CS(f(p,h ), O. 
and since 

/(/ (p . /'Mi)=/(/(/>,/i)./'), 


it follows that 
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This means that the point f(p,t x ) is recurrent as 
We thus see that if p is recurrent as /-*«, then all points 
f(p,t) will be recurrent as t-y ». Therefore it is natural 
to call the motion of f(p,l) recurrent as /-> ». The same 
remarks are also valid for points that are recurrent as 
t-y— «>. 

8 . Consider the set P of all recurrent points in the phase 
space of a mechanical system. For the following, it will 
be important to establish its measurability. We shall 
also need to consider the set P+{P~) of all points of phase 
space that are recurrent as ( 0 r t-y— «). 

Theorem 4.6. P + , P~, and P are sets of type G { and 
hence are measurable. 

Proof: Let 

Ti, r 2 , ..., r„, ... 

bean unbounded monotonic increasing sequence of positive 
numbers and 

??li • t ^7n> ■ • • 

a monotonic decreasing sequence of positive numbers 
which tends to zero as a limit. Let Pi denote the set of 
points p for which one can find a r>r„ such that 
p(/>, f{p^))<r) n . Consider a point p'ZPi, and let 
p(/ >/ */(/ >/ » T )) == <yj n - We then choose a positive number 
€ <(Vn —0)/2 and a $>0 sufficiently small so that for 
any point q for which p(p\q)<6, the inequality 
P{f(p f J)J(qS))< c is satisfied for all positive t^r. This 
is possible because of the continuous dependence of the 
phase motion on the initial point. We then have 
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p (?./>')<* < * < V 2 ( 1 i« — **)• 

P (/>',/(/>.*)) = ?. 

p(f(P',')> / (<7. *'))<«< Va ( t 4- — P>* 

and hence by the Triangle Property, it follows that 

P (?./(?. *))< V 

This means that from the condition p'£Pf t it follows that 

woe#, 

i.e., Pi is an open set. Obviously 

n p:, 

nml 

and this implies that P + is a set of type G 6 . 

The fact that P~ is also a set of type G t is proved 
analogously. 

The validity of the theorem for the set P then follows 
from the obvious relation 

p=p+hp-. 

Corollary. The sets E? —P and E —P are measurable . 

9. The theorems presented above were of a preliminary 
nature. We now proceed to prove the main theorems. 

Theorem 4.7. Let 9TI be a measurable set not containing 
points that are recurrent as >—»), where for each 
f(p,t) is defined for all t^Q (or t^ 0). If meas 9fR > 0, then 
one can find a bounded set 9Tl*C9tf of positive measure which 
is not recurrent as /->« (or ® ). 

We shall only prove the theorem for . The case 
tr ->—® is treated analogously. 
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Let 


Tli T3, ••• |T| | ,, , 

be an unbounded monotonic increasing sequence of pos¬ 
itive numbers, and 


TJh ^? 2 » • • • , Vity • • • 

a monotonic decreasing sequence of positive numbers 
tending to the limit zero. We let 317„ denote the points 
pdffZ such that 

f(P, t)ZS(p, Vh ) 

for all Since the set 5E does not contain any 

recurrent points, each point belongs to at least one 
of the 9H„ and, therefore, 


(4.i4) 

Let us show that 9TC* is measurable. To do this, we shall 
prove that each of the sets 9TC* is closed with respect to 911- 
Let p* be a point of 3TZ_which is a limit point for 9(11*. 
Let us assume that p*€ ( 3\l n ; then we can find a time 
t\>r n such that 

pO>*. f(P m . *i)) = P<u,. 

Now let € < (r) n —/9)/2, and let us chooseS >0 sufficiently 
small so that for any point q for which p(P*,q)< 5, the 
inequality 

P(/0>*, t(g*, *,))<«, 
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is satisfied. This is possible because of the continuous 
dependence of the phase motion on the initial point. 
Moreover, the point q may be chosen in SHI* since p* is a 
limit point of 3Tl„- Therefore, 

P (q, P*)<°< e <y O')* -- ?)* 


p(/>*» /(P*. 0) = ?* 

P *i)> /(?* 0)^* < ^'2’( 1 li» (*)• 

or by the Triangle Property 

Pfa. /(?. «i))<\. 

This contradicts the fact that q is in 911*. Therefore, 
p*£ 911* and 911* is closed. 

From (4.14), it follows that 

00 

S meas 9K„ ^ meas SD1 > 0, 

1 

and* therefore, we can find a set 9TC* such that 
meas 

Let x be a point of metric density of the set 911*, and 
consider the set 


ax;=-a«*n<s(*, y T *)- 

It is clear that 911* is measurable and that meas SH£*>0. 
The diameter of 911* does not exceed rj k /2 . But for our 
point p of 911*, we have 
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P(P<f(P, t)»Vk 
for alH > r k . Therefore, 

«t;n fW' k , 0=0 

for i>r h . Setting 9Tl* = 3TTk completes the proof of the 
theorem. 

Theorem 4.8. If A is a measurable set in E v such that 
any subset of it of positive measure is recurrent as £->oo 
(or £-> — 00 ), then all points of A, with the possible exception 
of a set of measure zero , are recurrent as t-+ qo (or oo). 

We prove the theorem by contradiction, and we assume 
that there exists a set STC of positive measure all of whose 
points are not recurrent as t-><» (or t->—a>). Then by 
Theorem 4.7, we can find a bounded measurable set 9TC* 
of positive measure which is not recurrent as 
(or t-^ —°o). However, this contradicts the conditions 
of the theorem. 

Let A be a measurable but otherwise arbitrary set of 
phase space of the dynamical system. We call the set A 
regular as t-+ oo if for any r >0, we can find a set A* such 
that meas A* < «» and 

U f(A,t)cA*. 

A set which is regular as t-> —® is defined in a similar 
way. 

Let B be a set and p a point of phase space; we call the 
point p recurrent with respect to the set B as (/->—«) 

if for any T> 0 , we can find a positive (or negative) 
number with \t f \>T, such that f(pZ) £B. 

Theorem 4.9. Let A be a regular set of points of phase 
space with invariant measure and W a set whose points are 
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recurrent with respect to A as t-y& (or t-y — 00 )• Then 
any measurable non-recurrent set 311* € PF w a set of measure 
zero. 

We prove the theorem by contradiction, and we assume 
that there exists a set 31 l*CW which is not recurrent and 
which has positive measure. Then we can give a positive 
number r such that the sets 

/(9R*. nr), n = 0, 1, ... 

are pairwise disjoint. Since A is regular, there exists a 
set A *, with meas A*< °° , such that 

U f(A, t)cA*. 

It is not difficult to see that for any we can 

choose a subsequence from the sequence f(p,nr ), O t l ,2, 

. . . , all of whose terms are points in A*. 

Let us consider the sets 

9K;=/(5r, nt), 

n as 0, 1, 2 . . . 

and let 

Each D n is measurable, being the intersection of two 
measurable sets. Moreover, the sets D n are disjoint. 
Therefore the series 

2 meas Z) v 0^-15) 


converges since 
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U D, C 4 *, 

v*0 

X 

meas U A. ^meas^ 4 * < cn 

v=0 

and 

® QO 

2 j meas/>v = measU A, <C 00 • (4 16 ) 

v=0 v * 

Further, 

/(£>„ — = — nx) ( 4 . 17 ) 

is the set of points p in SHI* for which we have 

f(p, nx)£A+. 

Each point of SHI* belongs to an infinite number of the sets 
(4.17) and, therefore, to each of the sets 

®' = & /(D «“ VX) - ( 4 . 18 ) 

M SB 1, I, ... 

Thus, for any n we have 

* 1 H* CH*. ( 4 . 19 ) 

By the invariance of measure, we obtain 

meas /(£>„ -vt) =meas Z) vn ( 4 . 20 ) 

and so by ( 4 . 18 ) and ( 4 . 20 ) 

m co 

meas 5 v„= 2 meas/ (^,* — vt ) = 2 meas D 

V3M vsl * 

Now, the series ( 4 . 15 ) converges and 
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2, meas D, 

M5St» 


is its remainder. Thus for any € >0, we can choose an 
integer N=N( e) sufficiently large that 

meas :1V w 2 (4.21) 


for all n'ZN. 

Taking into account relation (4.19), and inequality 
(4.21), we conclude that 

meas 3TI* < 3, 


and since e is an arbitrary positive number, this implies 
that 


meas flTC* =1). 


We have arrived at a contradiction, and this proves the 
theorem. 

10 . We now study the motion of a system of n grav¬ 
itating bodies in phase space making use of the analytical 
relations of celestial mechanics and the above general 
measure theoretical results on dynamical systems. 

Let us first consider the set Q of all points />€ H for 
which 


max {/*,,«} 00 as t —co. 

V 

max {r<,} -f* -c as / -► 4-00. 

•7 

We are not going to investigate the measure of Q. For our 
purposes, it will suffice to prove that Q satisfies the 
following purely logical theorem. 
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Theorem 4.10. From the assumption that 
meas Q > 0. 

it follows that all the points of Q, with the possible exception 
of a set of measure zero , satisfy the condition 

inf (min 

oiy 

Proof: Let Q~ denote the set of all points p€E Gn for 
which 


r* (/) =- max {r 4 j} oo as t —oo, (4.22) 


and let us show that this set is measurable. From condi¬ 
tion (4.22) it follows that for any point p€Q~ and any 
number 2?>0, we can find a positive number T so large 
that 

r* (t)^R for all t <— T. 


Let 


/*!, ^ 2 ’ 
p T 2 , 



be two unbounded monotonic increasing sequences of 
positive numbers. We let Q~{Rki Tm) denote the set of 
all points for which 

r* {t) ^ Rk for all t<^T m 

We now show that this set is closed. 

Consider the set G k of all points p€ E 6 * for which 

r * (P) R*' 

since r*(p) is a continuous function of the phase point, 
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Gk is an open set. Now arguing by contradiction, we 
assume that Q~{Rk,T m } is not closed. Then we can find 
a point p*€Q~{R k ,l m ) which is a limit point of this set. 
For this point, we can find a moment of time t'< — T m 
such that 

/(/>*, 06 G* 

Since the set G k is open, it Is possible to find a number 
e > 0 sufficiently small that 

0> s)CG t . 

Due to the continuous dependence of the phase motion on 
the initial point, we can choose a 8 > 0 sufficiently small 
that 


S), t')CZS(f(p*, /’), e)c6’*. (4.23) 


However, p* is a limit point of the set 
therefore we can find a point p' such that 


Q {RkTm} and 


T m ), 

P’ZSiP*, £). 


(4.24) 


From (4.23) and (4.24), it follows that 

/(/>'. 06 0 k , 

which leads to a contradiction. Thus the set Q~{R k ,T m } 
is closed. 

Let us next consider the set 


or=jji (<?-{*„ T m ) n S). (4.25) 

Its measurability follows from the results obtained above. 
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Now, Qr is the set of points of phase space each of which 
eventually leaves the set 

r *<-Kk (4.26) 

permanently as /-->—oo. 

From the obvious formula 


00 



it follows that the set Q~ is measurable. The set Q~ con¬ 
tains the points of phase space that each eventually leaves 
any set of the form 

r*(p)<R„ (4.27) 

fiasl, S, ... ' 

permanently as t —>— oo . 

Carrying out a similar argument for we can 

construct measurable sets 

Qt— U (Q+(R*, UflS). 

m=l 

analogous to the sets (4.25). Then 

is a measurable set consisting of the points in phase space 
each of which eventually leaves any set of the form (4.27) 
permanently as —« and which is also recurrent with 

respect to the set (4.26) as 
It is further clear that 

<?= u 

*=l 
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from which it follows first of all that Q is measurable. 

We are assuming meas Q> 0, and therefore, we can find 
a v such that for the set 

a=e-n(E-& + ) (4.28) 

we have 

meas Q v > 0. 


It is very important to note for the following that Q v is 
the set of points each of which, as t-+ —oo, eventually 
leaves any set of the form 

r* (P)$R n , 

n - 1 , 2 , ... 

permanently, and as t ~>» is recurrent with respect to 
the set 

Let 

an a.2» • • • * a*i, • • • 

be a monotonic sequence of positive numbers which con¬ 
verges to zero, and let A» denote the set of points p in 
E\ for which 

r * ( P ) = min {rij} ^ a n for all t ^ 0. 

Let us show that each A n is closed. Consider the set G * of 
all points for which 

r* (P)<a k ; 

since r#(p) is a continuous function of the phase point r 
G' k is open. Let us assume that A* is not closed. Then 
we can find a point />* £ A* which is also a limit point of 
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this set. For this point, we can then find a time /"> 0 such 
that 


e G' k . 

Since G' k is an open set, we can find an e>0 sufficiently 
small such that 


S(f(P*,t"))CG k . 

By the continuous dependence of the phase motion on the 
initial point we can choose a S > 0 so small that 

f(S(p*, 5), I") CS(M, n, 0 C G'k. (4.29) 

Now, p* is a limit point of A*. Therefore, we can find 
a point such that 


P" 6 A*, 
p"£S(p* } 5), 

From (4.29) and (4.30) it follows that 

f(p'\n£G' kt 


(4.30) 


and this leads to a contradiction. The set A* is closed 
and hence measurable. 

We now consider the proof of the theorem itself. 
We argue by contradiction. If the theorem is not 
true, then we can find a fx such that for the measurable set 


<?? = & A m , 


we have 


meas(#>0. 


(4.31) 
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Let us take an unbounded monotonic sequence of positive 
numbers 

• • • f Hm • • ■» 

and let us denote by W k the set of all points whose phase 
coordinates satisfy the inequality 


\H\ = 


y m ‘ ( x < ^ ~ v 



Since the left-hand side of this inequality is a function of 
the phase and continuous on H, each of the W h is an open 
set. Writing Q; as the sum of a denumerable number of 
measurable sets 


«=,u(«n»',) 

and taking into account inequality (4.31), we conclude 
that we can find a positive integer h such that for the set 

w=Qtnw h , 

we have 


meas W >0. 

For the following, it will be important to enumerate 
the properties of the points p in the set W: 

(1) if pZW, then as —co , /(/>,/) leaves any set of 

the form 

(?)<R . 

*1 = 1 , 2 , . . . 

permanently, i.e., ® as /->—oo ; 

(2) every point p £ W is recurrent with respect to the set 

r* (PXR, 


as /-> oo ; 
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(3) if p£W, then the motion of /(£,/) is such that 
r*(/)^ otp for all 0; 

(4) if p£W, then the motion of f(p,t) satisfies the 
inequality 


i*i= 


i=z\ 


U 


://i, 


for all real values of /. 

We let B denote the set of points p€ Jg for which 

min for all 0. (4.32) 

|Jf| < H h ; (4.33) 


From the preceding considerations, it follows that this set 
is measurable. Moreover, it is clear that if p€W } then 


f(p , 0 € j8 for all ^0 (4.34) 


so that B is invariant as $-*■«>. Let us show that the 
phase velocity v t considered as a function of the points of 
phase space, is bounded on the set B . The square of this 
velocity is given by 



(#)%(£)V(£)'+(£t+(£y+( 





(4.35) 


From inequality (4.32), it follows that the potential func* 
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tion U is bounded from above; therefore we can find an 
H*> 0 such that U^H*, and hence we have 

j 2 m i ( x t + y? + z i~) <.h-t h* , 

1=1 

i.e., the kinetic energy is bounded from above. Thus the 
functions #i 2 , y* 2 , and s',- 2 are bounded from above. Dif¬ 
ferentiating U with respect to x { , we obtain 


Since 


dU_ 

dn 



mimj 



dru 



we have that 



From this inequality and condition (4.32), it follows that 



is bounded from above. In a similar way, we can show that 
the functions 



are bounded. From these results and equation (4.35)» it 
follows that v\ and therefore also v , is bounded on the set 
Let 


A = Bf\S (n/? v ), 
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where S(nRy) is an open sphere with center at the origin, n 
being a positive integer (we are considering the motion of 
the n bodies in a coordinate system with origin at the 
center of mass). We have shown that the phase velocity 
is bounded on the set B . Therefore, we can find a positive 
number b such that |»| <b . Since A £ B } this estimate is 
also valid for the set A . Let us choose a r > 0 but other¬ 
wise arbitrary; for any p€B, the length of the path tra¬ 
versed by the point) f(p,t) |in phase space in the interval of 
time from 0 to r does not exceed br. Thus every point 
p£A is such that for all t in the interval 0 ^ r,we have 

f(p,i)€S(nR,+br). 

It then follows that 

U f(A, t)czS(nR, + fo) y 

and since S(nRy+br) is a bounded set, meas S(nR v +br) 
< oo and therefore, A is regular as . 

To complete the proof of the theorem, we must obtain 
a contradiction. 

Let us first show that there exists a set 911* C.W which 
is not recurrent as t ->—oo such that meas 911* > 0. Assume 
that this is not true. Then every measurable subset of 
W of positive measure is recurrent. However, in this case, 
we have by Theorem 4.8 that all points of W, with the 
possible exception of a subset of measure zero, are recurrent 
as / > oo. For the motion of f(p,t ), determined by any 
point p€W, we have that 

max {r^} -► oo as t —oo. 

This is impossible for recurrent points, and we arrive at a 
contradiction. Thus, there exists a set 9fH*CW^» with 
meas 911*>0, which is not recurrent as /-*—oo. By 
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Theorem 4.3, 9TC* is also not recurrent as 00 - Now, 
we know that each point p€W \s recurrent as /-*- 00 with 
respect to the set 

r*(p)<R^ 

Therefore, if at some moment of time the point /(/>,/') 
belongs to the set r*(p)<R yt the distances between the 
gravitating bodies will not exceed R v . Therefore, the 
distance of any of the bodies from the origin, located at 
the center of mass, does not exceed nR v , i.e.,/(£,0 £ S(nR p ). 
Hence, every point p that is recurrent with respect to the 
set r*(p)<R „, will be recurrent with respect to S(nR p ). 
These remarks, in conjunction with condition (4.34), allow 
us to say that every point p£ W will be recurrent as <» 
with respect to the set 

A = linS(nIiJ, 

a set that is regular as /-><». We have thus arrived at 
the conditions of Theorem 4.9 from which it follows that 
meas 9E* = 0.|Our contradiction is obtained, and the 
theorem is proved. 

We shall call a system of n gravitating bodies almost 
stable as oo if 

max {r,,} -I—*oo as t 4-oo, 
iJ 

and if we can find a number a > 0 such that the condition 

min {rtf} > a 

iJ 

is satisfied for all t^ 0. 

We shall call a system of n gravitating bodies stable 
as tr+co if we can find positive numbers R and a such that 
for all t^O 

max {tv} < H, min {/■„} > a 
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At present, it is not known whether there exist almost 
stable systems which at the same time are not stable. 

We shall say that a system of n gravitating bodies 
contains a body which has come from infinity if 

max {r^} co as* t -► —oo. 

v 

Theorem 4.11. For all initial conditions , with the pos¬ 
sible exception of a set of initial conditions of measure zero , 
a system of gravitating bodies which contains a body that has 
come from infinity cannot be almost stable . 

We prove the theorem by contradiction. 

If the theorem is not true, we can find a number a>0 
such that the set Q* of all points p satisfying the 
conditions 

max {r >v } -> oo as t —oo, 

7 

max {r lV }-► oo as 

i; 

inf(min{r„})>a 

ij (4.36) 

is a set of positive measure: 

meas Q*^>Q. (4.37) 

However, Q*C.Q » where Q is the set of points E 
satisfying conditions (4.35). From the fact that Q*€Q 
and the condition (4.37), it follows that 

meas Q >0. (4.38) 

Now, by Theorem 4.10, we know that when condition 
(4.38) is satisfied by all points of Q , with the possible 
exception of a subset of measure zero, the condition 



96 QUALITATIVE METHODS IN THE MANY BODY PROBLEM 


inf (max (r iy )) = 0. 

0 ij 

must be satisfied. 

This result, clearly, is also valid for the subset Q *, 
and this contradicts condition (4.36). 

The theorem is thus proved. 



CHAPTER 5 


Analysis of Some Cases of the Evolution 
of a System of Gravitating Bodies 

1. The two-century-old development of the planetary 
cosmogony has often been described as a haphazard 
alternation of diverse hypotheses, as a process of trial and 
error that has produced no definite results. This, how¬ 
ever, is not so. The development of cosmogony as a 
science is a process that has followed the general pattern 
of growth of knowledge. Cosmogonical knowledge has 
been acquired in the form of relative truths, and as these 
truths have evolved, we observe the gradual formation 
of the elements of an objective science, which thus be¬ 
comes continually more perfect. Of course, the growth 
of cosmogony, like that of other sciences, has been a very 
complicated and contradictory thing, in which true and 
false principles have emerged side by side. However, 
in the course of its development, certain ideas have 
appeared that have persevered, and with modification and 
expansion, have remained as part of the science in one form 
or another. The history of science shows that such evolv¬ 
ing ideas constitute that part of the science which, to 
the extent of its growth, form the elements of objective 
knowledge. 

In the cosmogony of the planets, one such idea is, 
first of all, that the solar system and naturally formed 
systems similar to it evolved out of rotating matter 
scattered in space. Much controversy arose about the 
physical state of the matter, about its genetical relation¬ 
ship to the various other cosmical bodies, and about the 

97 
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nature of the process by which the matter was accumulated 
and transformed into large bodies. However, a general 
basis has stood the test of time and is being confirmed 
by the experience of the two-century-old development of 
the planetary cosmogony. This basis thus comprises the 
elements of objective knowledge. The idea that the 
planets were formed from rotating matter scattered in 
space follows almost inevitably from the properties of 
the solar system and the natural laws governing it. 

Thus the basic questions we have to consider in cos- 
mogonical investigations are the process by which the 
matter scattered throughout wide regions of space was 
accumulated in a small volume, the process by which this 
matter was condensed, and the process by which tiny 
particles were unified into large bodies when the cosmical 
systems were formed. There is no doubt that gravita¬ 
tional forces play an essential role in these phenomena. 
However, today, it is regarded as well established that 
gravitational forces alone are insufficient to explain their 
occurrence, and of fundamental importance also are phys¬ 
ical phenomena involving energy conversion. Therefore, 
the full analysis of these cosmogonical processes is beyond 
the scope of gravitational mechanics. 

Nevertheless, the theorems of celestial mechanics can 
be of considerable help in these investigations. We begin 
by considering the following. The rotating matter scat¬ 
tered in space, from which the solar system was formed 
and which we take as its initial state, i.e., as the object 
whose evolution is impending, can be regarded as a system 
of many gravitating bodies. The final product of evolu¬ 
tion, the solar system, is also a system of many gravitating 
bodies which are fewer in number and larger in size. 
Without going deeply into the processes and mechanism 
by which the initial system was transformed into the final 



EVOLUTION OF A SYSTEM OF GRAVITATING BODIES 99 


one, but by comparing them only as two gravitating 
systems of bodies, we shall discover that this transforma¬ 
tion cannot be realized solely on the basis of the laws of 
gravitational mechanics. By looking only at the mechan¬ 
ical side of the question, we cannot, to be sure, explain the 
physical factors which are important to the transformation 
under consideration. Nevertheless, a consideration of the 
mechanical side of the question does yield significant 
results. We may not be able to discover the nature of 
these physical factors, but we can get some idea about 
them by determining the mechanical consequences of these 
factors. 

2. Let us consider n gravitating bodies P 0 , Pi, , 
P*_i, and let denote the distance between Pi and Pj. 
We shall say that a completely hyperbolic approach of the 
bodies occurs if all the distances satisfy the condition 

00 as /->— oo 

It is not difficult to see that a completely hyperbolic 
approach is a type of motion that is actually possible for a 
system of n gravitating bodies. Namely, we know that 
in purely dynamical problems, it is legitimate to reverse 
the sign of the time. Therefore, if the inequalities (3.24) 
hold for the case of time varying negatively, then from the 
appropriate reformulation of Theorem 3.4, it will follow 
that a completely hyperbolic approach takes place. More¬ 
over, since r u and r\j are continuous functions of the 
coordinates of phase space of the system of gravitating 
bodies, the set of points of phase space whose coordinates 
satisfy condition (3.24) is an open set. Hence, it follows 
that the set of points of phase space which are initial 
points for the “regime” of a completely hyperbolic ap¬ 
proach is a set of positive measure. 
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Let us now show that the motion in which the cond i tions 
fiji<x> as /-> — oo (5-1) 

and 

rij(t)<R as 0, (5.2) 

are satisfied simultaneously, where R is any suitable 
positive number, is impossible for a system of n gravitating 
bodies. 

Disregarding a set of initial states of measure zero (the 
case H = 0), we can assert that from (5.1) and Theorem 1.1, 
it follows that H> 0. Therefore, by Theorem 1.2 at least 
one of the distances between the bodies tends to infinity 
as , which contradicts condition (5.2). 

We thus see that for a completely hyperbolic approach 
of gravitating bodies, a system cannot arise in which the 
mutual distances remain bounded for all 0. However, 
if the bodies approach each other, then the unification of 
not all, but some portion of them into a subsystem with 
bounded is still possible. This was rigorously proved 
for the three-body problem in the theory of capture de¬ 
veloped by 0. Yu. Schmidt, G. F. Khilmi, N. N. Parisky, 
and G. A. Merman. For a larger number of bodies, it 
should be possible to give examples in which such motions 
occur, however, they would be somewhat artificial and 
would entail very formidable computations. At any rate, 
it is possible to get a sufficiently clear picture of our 
phenomena by means of the following considerations. 

If from a number of bodies that have come from 
infinity, a stable subsystem is isolated following the ap~ 
proach of the bodies, then on the one hand, J7->0 as 
/->— co, and on the other hand, there are distances be¬ 
tween the bodies that are bounded on the set t ^ 0. There- 



EVOLUTION OF A SYSTEM OF GRAVITATING BODIES 101 


fore, we can find a constant U 0 and a time T such that 
U> U 0 for t> T. 

Then from the formula 

t 4 

I 1 (0 = II + II' t -I- 2M J fu+2H) dtzdh, 
rr ~v~ 

in which we must take JT> 0 , we see that following the 
approach of the bodies, the moment of inertia P(t) will 
increase more rapidly as than as —» despite the 
fact that all the distances increase indefinitely as /->—<», 
and only a portion of them increase indefinitely as /->». 
This means that the gravitational association of a portion 
of the bodies is accompanied by the rise of motions 
involving increased velocities. 

The physical interpretation of these facts is as follows. 
The too high kinetic energy which the bodies have as the 
result of approaching from infinity prevents their complete 
gravitational association. However, when certain condi¬ 
tions are satisfied, a portion of the bodies can form a stable 
subsystem. This can happen when the gravitational in¬ 
teraction of the bodies as they approach each other results 
in a redistribution of the energy causing some parts of 
the system to begin to move rapidly with respect to the 
others. If the kinetic energy of the relative motion of 
parts of the system constitutes a large percentage of the 
total energy, these parts can then become stable systems. 
The unification of a portion of the bodies into one or 
several stable subsystems is necessarily accompanied by 
the “draining off” of the energy in the system, i.e., by 
the appearance of relative motions of the subsystems in 
which the surplus of energy preventing any association of 
the bodies is spent. Theorem 3.8 shows that if at least one 
of the bodies in a system of n gravitating bodies has 
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sufficiently high kinetic energy, the other bodies cannot 
comprise a dissipative system. 

We thus see that the association and scattering of 
gravitating systems are connected in the closest way. 
If one of these phenomena occurs, it will be compensated 
for by the occurrence of the other. Association and 
scattering in gravitating systems are thus divergent aspects 
of the one phenomenon of the gravitational interaction of 
matter . 

The above results do not mean that whenever we have 
a hyperbolic approach of many bodies, it is impossible 
for all of them to be unified into a stable system under 
their mutual attraction. Our investigations merely point 
out that this is impossible as a purely dynamical process. 
However, if we make allowances for processes of a more 
general type, then this may conceivably happen. 

In fact, if whenever the bodies drawing near under 
their mutual attraction gives rise to or is accompanied by 
some physical process in which the surplus of kinetic 
energy preventing their association is absorbed, then 
association can be realized without the phenomenon of 
scattering or else with a relaxation of it. In this case, 
scattering is “replaced” by a physical process that absorbs 
the surplus of mechanical energy and converts it into 
non-mechanical form. 

For meteoric dust-clouds, the most important form of 
“drainage” of mechanical energy occurs when a portion 
of it is converted into thermal energy during non-elastic 
collisionsof the particles. When there is such a “drainage,” 
association is conceivable without the scattering of dust 
particles or with the relaxation of scattering. In par¬ 
ticular, it is possible for meteoric dust particles to be 
captured when a star passes a meteoric dust-cloud closely. 
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The dynamics of stellar passages through a cloud of 
meteoric dust particles was investigated by T. A. Agekyan 
who worked out a theory of capture tied up to the particles’ 
loss of kinetic energy in mutual nonelastic collisions. 

3. We can formulate some further general results by 
using Theorem 4.11 which was proved in Chapter 4. To 
do this, we must know under what circumstances the 
evolution of any particular system of real bodies is 
describable by means of the dynamical laws for a system of 
n gravitating bodies. 

Real bodies such as the stars, planets, small solid 
particles, etc., are idealized as mass points. However, the 
idealization of real bodies as mass points is permissible 
as long as the distances between them are sufficiently great 
in comparison to their dimensions. Once the bodies are 
distributed so closely that the distances between them are 
comparable to their dimensions this idealization ceases 
to be correct. In this case, various physical phenomena 
arise which are accompanied by the conversion of mechan¬ 
ical energy into a non-mechanical form. Finally, when 
approaching closely, these bodies may undergo considerable 
changes: they may collide or be disintegrated or, con¬ 
versely, they may be unified, etc. Therefore, whenever 
the motion of real bodies is to be described in a given 
interval of time by the dynamical equations for a system 
of gravitating bodies, we must assume that the distances 
between the bodies do not become smaller than some 
positive number a whose value will depend on the concrete 
properties of the given system of bodies. 

With these remarks, let us consider a system of grav¬ 
itating bodies Po, Pi, ... , P„—i, which we assume to have 
the following properties: 
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r 0j (0 00 as 00 » 

y=i.2,.. 

r <y(04"*°° M l-*4-°o t 

5=1,2.n-1 

i.e., we consider a system of bodies Pi, Pf, • • • , P*-i 
which does not get scattered as /-><», to which is adjoined 
the body Po coming in from infinity. We next show that 
if a physical system is formed after Po is adjoined to the 
system consisting of Pi, Pi, ... , P*_i which is stable as 
/->oo, then the process by which this system is formed 
could not be purely mechanical: it is necessarily accom¬ 
panied by the conversion of mechanical energy into non- 
mechanical form. 

In fact, Theorem 4.11 implies that except for a set of 
initial conditions of measure zero, for every system of 
many bodies containing a body incoming from infinity 
but none outgoing to infinity, the greatest lower bound of 
the mutual distances r i3 is zero. Therefore, arbitrarily 
close approaches of the bodies in the system will occur, 
which means also the nonmechanical phenomena that 
accompany them. 

For the sake of simplicity, we have considered the case 
where the system contains only one body coming in from 
infinity. However, our conclusions still hold for any 
number of such bodies. 

4# We now derive an inequality which we shall use 
below. 

We consider n gravitating bodies Pi, Pi,...»P« with 
masses mu m i, ...»m#, and we describe their motion 
with respect to a set of coordinates with origin at the 
center of mass of the system. 

Consider the following identities whose validity can 
easily be verified directly: 
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=[M - v/.'f+M - MflP +M - *X,n + rf. 

1 » 2 , ...» x 

Multiplying these identities respectively by mi, m 2 , . . . , 
and w n and summing, we obtain 


2 T=P+U\ (5.3) 

where T is the kinetic energy and P and Q 2 are given by 

H 

7>= IS — y ^y + 


+2 ‘t _ z ^) 2 +IStt - w 

<=1 * tel < 


(5.4) 



(5.5) 


Now, consider the first sum in the expression on the 
right-hand side of (5.4). The variables x if y i} x f i} and y\ 
appearing in this sum are not independent. They are con¬ 
nected by relations (1.6) expressing the area integrals. If 
we use Lagrange multipliers to determine the relative mini¬ 
mum of the first sum in equation (5.4), we find that it is 
equal to c\/J 2 . 

In an analogous way, we find that the relative minima 
of the two other sums on the right-hand side of (5.4) are 
respectively 


c\ and c\ 
jt jr 
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Using these results and recalling that 
Q 2 — c 3> 

we find that 


P> 


©i 
J 2 ' 


(5.6) 


From relations (5.3) and (5.6) it follows that 


2 T> + 


or 


j 2 > 


e 2 

2 t-w ' 


From the energy integral, we have that 


(5.7) 


2T = 2{U+H)\ 


substituting this expression for IT in (5.7), we obtain the 
required inequality 




e 2 

2(U+H)-W 


(5.8) 


5. A system of a large number of gravitating bodies is 
called stationary, if despite the movement of its component 
bodies, the distribution function for the distances between 
the bodies remains practically constant and in this sense 
is an invariant of the motion. Clearly, a system can be 
stationary only if its total energy is negative. 

If a system is stationary, then 


dt * 


/•'=const; 


0. 
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Therefore, it follows from the Lagrange-Jacobi equation 
that 

£/+2tf=0. 

From this equation and the energy integral, we have that 

T=-H, U = —2//, 

or taking into account that if< 0 , we can write 

T=\H\ } U—2\H\. (5.9) 

6 . Let a physical system 5 be in the state (A) at time 
to in which it is a collection of n 0 mutually gravitating 
bodies. Suppose that at time t\>k , the system 5 passes 
into a new state (B) in which it is a system of wi gravitating 
bodies. As to the transition of S from the state (4) to 
the state (B), we assume nothing other than it has occurred 
as the result of the gravitational interaction of the moving 
bodies. Of course, of essential importance in this transi¬ 
tion are also the nonmechanical phenomena which are not 
manifested in the dynamical laws for gravitating bodies. 
Under this formulation of the problem, the number of 
bodies n 0 and wi in the respective states (^4) and ( B ) are 
not necessarily equal. 

We let 

7o, u 0l i//oi, ji eg, Qg 

and 

T u U U \H U J\,& U Of 


denote the values of the quantities 7\£7, |//|, 7 2 , fi 2 , and 
0 2 for the respective states (^4) and (5) (in other words, 
at the times to and * 1 ). 

Suppose that state C<4) satisfies the condition 
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d*JX 



(5.10) 


and suppose that after transition into state (. B), S is a 
stationary system of bodies occupying a region so small 
that 




A* 

w o 


2 1// 0 | — aj 


(5.11) 


Let us study the transition of S from state (A ) to state 
(B). From the Lagrange-Jacobi equation (1.14) and con¬ 
dition (5.10), it follows that 

f/„ + 2// u ^0. 


From this inequality and the energy integral, we obtain 
Taking into account that 7/<0, we have that 


\H„ \>T„ 

or that 

2|7/ 0 | > 7 (> . 

Subtracting S2§ from the left and right-hand sides of this 
inequality, we obtain 

2 1// 0 |—2g>7’ 0 —ar; 

Now, from the definition of 7» and S2§, it follows that 

o. 


and therefore that 
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Thus if 0§>O, the following inequality is satisfied: 






3>0. 


(5.12) 


This inequality shows that condition (5.10) is math¬ 
ematically consistent (admissible). 

With the system 5 in state ( B ), we now write inequality 
(5.8) using the condition 


V , = 2|/f,l. 


which holds because of the stationarity of state (B). We 
obtain 


/?>: 


21//, | — tif 


2 ' 


(5.13) 


From inequalities (5.11) and (5.13), it follows that 


<r,^ 


H 0 


2|/f,l-Uf 2|// 0 |-12 S 


Now, in order for this inequality to hold, at least one of 
the following inequalities must be satisfied: 



(5.14) 

e*<eS, 

(5.15) 


(5.16) 


The results obtained, of course, do not permit us to 
clarify the character of the non-mechanical processes that 
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occur in the transition of system 5 from state (-4) to state 
(j B). However, some information concerning these proc¬ 
esses can, nevertheless, be obtained. Namely, it is pos¬ 
sible to indicate some of the mechanical consequences of 
the physical phenomena that accompany the transition 
of S from state (. A ) to state (. B ). 

Thus, let us consider inequalities (5.14), (5.15), and 
(5.16). If inequality (5.14) is satisfied when S passes 
from state (A) into state ( B ), this means that the transi¬ 
tion is accompanied by phenomena that cause a decrease 
in the mechanical energy of the system. Thus, for ex¬ 
ample, if the system in the state (.4) is a cloud of solid 
particles, then the nonelastic collisions between these 
particles will cause a decrease in the mechanical energy 
of the system since as the result of these collisions a portion 
of the kinetic energy is converted into heat. 

If, however, inequality (5.15) is satisfied, this means 
that the transition from state (4) into state (. B ) is accom¬ 
panied by phenomena for which the moment of orbital 
momentum of the bodies comprising the system S in the 
state (A) passes into the moment of axial rotation of the 
bodies which make up 5 in the state ( B ). The transition 
of the orbital moment into the axial moment is possible 
if there are collisions between the bodies, or if many 
bodies are unified into a single larger body. 

We must still consider the case where inequality (5.16) 
is satisfied. If all the orbits of the bodies in the system 
were circular, fl 2 = 0. As the orbits deviate more from 
circles, i.e., as the radial velocities r' become larger, the 
value of ft 2 will become larger. 

Therefore, we can regard Q 2 as a measure of the 
deviation of the actual orbits from circular orbits. But 
then the fact that (5.16) holds, means that the transition 
of 5 from state (.4) to state ( B ) is accompanied by a 
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rounding off of the orbits, 1.0., by the orbits changing into 
a more circular shape. This can occur as the result of the 
interactions and disintegrations of the bodies when they 
collide, or when many bodies are unified into a single 
larger body. 

It is understood, of course, that one cannot exclude the 
possibility that the phenomena may be so related that 
after the transition of the system 5 from state {A) to state 
(-8), all three inequalities (5.14), (5.15), and (5.16) will be 
satisfied. 

The analysis of the problem in question by means of 
the methods of celestial mechanics does not permit us to 
obtain more concrete results. To this end, it is necessary 
to carry out further cosmogonical investigations. How¬ 
ever, we must point out that good agreement does exist 
between our results and those given by present-day 
cosmogony. The investigations of academician 0. Yu. 
Schmidt, confirmed later on by the work of L. E. Gurevich, 
A. I. Lebedinskii and other research workers, have shown 
that the formation of the planetary system from scattered 
matter rotating about the sun was accompanied by con¬ 
versions of the kinetic energy of this matter into thermal 
energy, as well as by the unification of the scattered matter 
into the planets, and by the rounding out of their orbits. 
The conversion of mechanical energy into nonmechanical 
form turned out to be the basic process, and the two other 
processes its consequences. 
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